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Using  9.0  fb~^  of  data  collected  with  the  CLEO  II.V  detector  at  the  CESR 
storage  ring,  we  compare  the  lifetimes  of  neutral  D mesons  decaying  via  — )•  K~K^ 

and  — >■  to  to  measure  the  mixing  parameter  ycp- 

We  find  ycp  — —1.1  ± 2.5  ± 1.6%,  where  the  first  error  is  statistical  and  the  second 
is  systematic.  This  corresponds  to  a 95%  confidence  interval  —6.7  < ycp  < 4.5%. 
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CHAPTER  1 
THEORY 


In  this  document,  we  describe  an  analysis  of  meson  lifetimes  when  the 
decays  into  CP'^  final  states.  In  this  first  chapter,  we  discuss  the  physics  behind  this 
work  and  describe  why  this  particular  analysis  is  interesting. 

1.1  Particle  Physics  and  the  Standard  Model 

Particle  physics  is  in  many  ways  the  quest  to  find  the  answers  to  just  two 
questions:  “What  are  the  fundamental  building  blocks  of  matter?”  and  “How  do  these 
building  blocks  interact  with  each  other?” . In  the  last  thirty  years  our  understanding 
of  the  universe  has  advanced  so  tremendously  that  we  now  have  partial  answers  to 
these  two  questions.  We  now  have  a convincing  picture  of  the  fundamental  structure 
of  all  observable  matter  in  terms  of  a small  set  of  building  blocks,  the  elementary 
particles.  The  extremely  successful  (if  incomplete)  theory  describing  the  behavior  of 
and  interactions  among  these  building  blocks  is  called  the  Standard  Model. 

First  we  describe  the  fundamental  constituents  of  matter  and  the  forces  they 
interact  through,  according  to  the  Standard  Model.  In  the  Standard  Model,  three 
generations  of  leptons  and  quarks  are  supposed  to  be  the  fundamental  building  blocks 
of  matter.  Each  generation  of  quarks  and  leptons  has  two  components.  For  the  quark 
generations,  shown  in  Table  1.1,  there  is  a quark  with  electric  charge  +2/3  (among 
other  quantum  numbers)  and  a partner  quark  whose  electric  charge  is  —1/3.  For  the 
lepton  generations,  shown  in  Table  1.2,  there  is  a lepton  with  electric  charge  —1  and 
a partner  lepton  (a  neutrino)  whose  electric  charge  is  0.  The  quarks  and  leptons  are 
fermions  (particles  with  half  odd  integer  spin;  natural  units  are  used  throughout  this 
document). 
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Table  1.1: 

Quark  generations  and  some 

of  their  properties. 

Name/Flavor 

Symbol 

Electric 

Charge 

Weak 

Isospin 

Has 

Color 

Mass  (MeV) 

Up 

u 

+2/3 

+1/2 

Yes 

2-8 

Down 

d 

-1/3 

-1/2 

Yes 

5-15 

Charm 

c 

+2/3 

+1/2 

Yes 

1000-1600 

Strange 

s 

-1/3 

-1/2 

Yes 

100-300 

Top 

t 

+2/3 

+1/2 

Yes 

?«175000 

Bottom 

b 

-1/3 

-1/2 

Yes 

4100-4500 

Table  1.2:  Lepton  generations  and  some  of  their  properties. 

Name 

Symbol 

Electric 

Charge 

Weak 

Isospin 

Has 

Color 

Mass  (MeV) 

Electron  Neutrino 

0 

+1/2 

No 

< 0.000015 

Electron 

e 

-1 

-1/2 

No 

0.511 

Muon  Neutrino 

0 

+1/2 

No 

< 0.19 

Muon 

-1 

-1/2 

No 

105.7 

Tau  Neutrino 

Ur 

0 

+1/2 

No 

< 18.2 

Tau 

T 

-1 

-1/2 

No 

1777.0 

The  quarks  and  leptons  make  up  all  matter;  each  quark  and  lepton  has  a corre- 
sponding antimatter  partner  that  has  opposite  quantum  numbers,  but  the  same  mass 
as  the  matter  particle.  Ordinary  matter  (atoms  and  molecules)  is  composed  of  only 
a subset  of  elementary  quarks  and  leptons.  Two  quarks  (up  and  down)  compose  pro- 
tons and  neutrons,  and  one  lepton  (the  electron)  makes  up  the  matter  that  humans 
encounter  in  everyday  life.  Except  for  their  masses,  each  quark  or  lepton  generation 
is  a duplicate  of  another. 

The  forces  between  the  elementary  particles  are  mediated  by  the  gauge  bosons, 
shown  in  Table  1.3,  and  all  have  spin  1.  The  Standard  Model  describes  three  kinds 
of  forces.  The  electromagnetic  force,  which  is  familiar  to  us  all,  is  mediated  by  the 
photon.  The  electromagnetic  force  is  responsible  for  binding  nuclei  and  electrons  into 
atoms,  for  example.  Only  particles  with  electric  charge  can  interact  electromagneti- 
cally.  The  strong  force,  mediated  by  the  gluon,  couples  particles  that  have  color,  and 
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Table  1.3:  Gauge  Bosons  and  some  of  their  properties. 


Name 

Symbol 

Force  Mediated 

Electric  Charge 

Mass  (GeV) 

Photon 

7 

Electromagnetism 

0 

0 

Gluon 

9 

Strong 

0 

0 

Weak 

0 

91.187 

w+- 

w+- 

Weak 

±1 

80.40 

is  responsible  for  holding  quarks  together  inside  hadrons.  Color  is  a kind  of  “strong 
charge,”  very  much  analogous  to  electric  charge.  Each  quark  flavor  comes  in  three 
different  colors,  and  a total  of  8 gluons  are  required  to  describe  the  strong  forces 
between  quarks.  Leptons  have  no  color,  and  therefore  do  not  interact  via  the  strong 
force.  This  is  one  of  the  primary  differences  between  quarks  and  leptons.  The  weak 
force  is  mediated  by  the  W'^,  W~,  and  Z^.  Particles  with  “weak  charge,”  or  “weak 
isospin,”  can  interact  via  the  weak  force.  An  example  of  a weak  interaction  is  the 
^-decay  of  heavy  nucleii. 

Note  that  in  Table  1.3  there  is  no  force  carrier  for  the  most  familiar  force, 
gravity.  This  is  because  gravity  is  not  included  in  the  Standard  Model,  and  is  indeed 
one  of  the  reasons  the  Standard  Model  is  incomplete.  Since  the  masses  of  the  particles 
are  so  small,  this  turns  out  to  be  not  important  in  particle  physics. 

Quantum  chromodynamics  (QCD)  is  the  part  of  the  Standard  Model  that  de- 
scribes the  strong  force;  i.e.  the  interactions  between  quarks  and  gluons.  It  explains 
why  free  quarks  have  never  been  observed,  by  requiring  that  quarks  be  confined  in 
“color  singlets,”  that  is,  objects  that  have  no  net  color.  The  possible  colors  for  quarks 
are  r (red),  b (blue),  and  g (green).  Anti-quarks  have  “opposite”  colors:  f (anti-red), 
b (anti-blue),  and  g (anti-green).  Colorless  objects  that  can  be  constructed  include 
mesons,  qq  states  where  the  quark  and  anti-qaurk  colors  cancel  each  other  (for  ex- 
ample, a red  quark  and  an  anti-red  anti-quark),  and  baryons,  which  are  three  quark 
states  where  each  quark  is  a different  color  (red  -I-  green  -f-  blue  = white  = “color- 
less”). Together,  mesons  and  baryons  (and  anti-baryons)  make  up  a class  of  particles 
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called  hadrons.  Quarks  and  anti-quarks  are  confined  into  hadrons  in  QCD  with  a 
steep  inter-quark  potential  that  becomes  extremely  strong  at  large  distance.  Using  a 
meson  as  an  example,  if  we  put  energy  into  the  system  to  try  to  separate  the  quark 
and  anti-quark,  at  some  point  it  becomes  energetically  favorable  to  produce  a new 
quark  pair.  Two  color  singlet  mesons  result,  and  we  never  actually  see  any  free  quarks. 
On  the  other  end  of  the  distance  scale  when  our  quark  and  anti-quark  are  close  in 
our  meson,  the  inter-quark  potential  is  smaller  and  the  quarks  are  comparatively  free 
to  bounce  around. 

The  Standard  Model  is  very  successful  in  explaining  all  existing  experimental 
data,  despite  the  lack  of  treatment  for  gravity.  However,  several  additional  issues 
remain  unclear.  For  instance,  the  Standard  Model  mechanism  by  which  particles 
acquire  mass  involves  the  existence  of  a yet  unobserved  boson,  the  “Higgs”  particle. 
The  Standard  Model  says  nothing  about  what  its  mass  may  be.  However,  experi- 
mental evidence  suggests  it  must  be  below  200-300  GeV.  Why  is  that?  Further,  the 
Standard  Model  also  provides  a mechanism,  CP  violation,  that  could  explain  the 
excess  of  matter  over  anti-matter  in  our  universe.  ^ However,  is  there  “enough”  CP 
violation  allowed  in  the  standard  model  to  allow  creatures  like  ourselves  to  exist,  or 
are  there  some  still  unknown  processes  required  for  that  to  happen?  There  are  many 
other  questions  about  which  the  Standard  Model  is  silent.  All  told,  there  are  some 
twenty  free  parameters  in  the  Standard  Model  (depending  how  you  count),  which 
could  imply  that  the  Standard  Model  is  a low-energy  approximation  of  a greater  the- 
ory that  can  answer  questions  the  Standard  Model  cannot.  Therefore,  experimenters 
make  measurements  that  can  be  checked  for  consistency  with  Standard  Model  pre- 


^Were  this  not  true,  none  of  us  would  be  here. 
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dictions.  An  observed  inconsistency  would  be  a signature  for  “New  Physics”  beyond 
the  Standard  Model,  which  could  result  in  answers  to  some  of  these  questions. 

1.2  Symmetry  Operations  C,  P,  and  CP 

Symmetry  operations  upon  a system  transformation  the  system  in  a way  which 
leaves  it  indistinguishable  from  the  original  configuration.  Such  symmetries  manifest 
themselves  in  the  laws  of  physics. 

Charge  conjugation  is  one  such  symmetry  operation.  The  charge  conjugation 
operator,  C,  is  a unitary  operator  {C~^  = C^)  changes  a particle  a into  its  anti-particle 
a and  vice-versa.  Therefore,  we  can  write 

C\a)  = 'qc\a)  (1.1) 

where  t]c  is  a phase  factor  Certainly,  we  can  let  C operate  on  |o)  twice  to  recover 
our  original  state  |a) 

C^|a)  = T]c{a)C\a)  = r]c{a)r]c{a)\a)  = |a)  (1.2) 

which  shows  that  rjc{a)T]c{a)  = 1,  and  it  follows  that  ric{a)  = r]c{a)*.  If  (and  only 
if)  a is  its  own  anti-particle,  then  it  is  an  eigenstate  of  C ^ and  T]c{a)  = rjc{a)  = ±1. 
Examples  of  such  eigenstates  are  mesons  composed  of  a quark  and  anti-quark  of  the 
same  flavor  (qq),  such  as  the  7t°,  (j),  and  p°.  The  eigenvalues  of  C for  qq  mesons  are 
qc  — (~1)^^'^)  where  L is  the  meson’s  orbital  angular  momentum  quantum  number 
and  S is  its  spin.  An  eigenstate  with  an  eigenvalue  of  -1-1  is  said  to  be  even  under 
charge  conjugation;  if  the  eigenvalue  is  -1,  it  is  said  to  be  odd.  Eigenvalues  of  the 


^If  an  equation  can  be  written  as  A|a)  = A|a),  the  state  |o)  is  said  to  be  an  eigenstate  of  the 
operator  A,  and  A is  its  eigenvalue. 
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Table  1.4:  Eigenvalues  of  C and  P for  some  qq  mesons. 


Meson 

c 

p 

7T^ 

+1 

-1 

V 

+1 

-1 

-1 

-1 

U 

-1 

-1 

<t> 

-1 

-1 

charge  conjugation  operator  for  some  of  the  mesons  used  in  this  analysis  are  listed  in 
Table  1.4. 

Another  example  of  a symmetry  operation  is  parity.  The  parity  operator,  P,  is 
a unitary  operator  that  reverses  the  sign  of  an  object’s  spatial  coordinates,  such  that 
P'tp{t,r)  = ip{t,  —f).  For  a particle  with  four  momentum  p = {E^p)  and  helicity  A, 
the  parity  transformation  can  be  written 


P|a(p,  A))  = T]p\a{p,  -A))  = r}p\ap) 


(1.3) 


where  p = (P,  —p)  and  |ap)  is  the  parity-transformed  version  of  |a).  Parity  reverses 
the  momentum  direction  of  |a),  but  not  its  spin;  therefore,  the  helicity  \ — s-p 
changes  sign  under  parity.  Similar  treatment  as  for  C shows  that  eigenvalues  of  P 
must  also  be  ±1.  Mesons  are  eigenstates  of  parity,  with  eigenvalues  rjp  = (— l)^'*'^, 
where  L is  again  the  meson’s  orbital  angular  momentum  quantum  number.  We  say 
an  eigenstate  has  even  parity  if  its  eigenvalue  is  +1  and  odd  parity  if  its  eigenvalue  is 
— 1.  Eigenvalues  of  the  parity  operator  for  some  of  the  mesons  used  in  this  analysis 
are  listed  in  Table  1.4. 

The  CP  operation  is  then  the  application  of  the  parity  and  charge  conjugation 
operations  to  a system;  that  is,  a particle  turns  into  its  anti-particle  while  at  the  same 
time  the  system  is  reflected  and  inverted.  Eigenstates  of  C and  P,  like  qq  mesons,  are 
also  eigenstates  of  CP,  with  eigenvalues  pcp  — VcflP  for  single  particle  systems;  for 
systems  of  particles,  we  pick  up  an  additional  multiplicative  factor  of  (— 1)^  where  I 
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is  the  orbital  angular  momentum  of  the  system.  For  future  reference,  the  neutral  D 
mesons  of  interest  in  this  document  transform  thus  under  C,  P,  and  CP 


C|po  > = 

|P0  > 

C|po  > = 

|P°  > 

P\D^  > = 

-IP*^  > 

P|P0>  = 

-|P0  > 

(1.4) 

cpipo  > = 

-|P0  > 

CP|po>  = 

- > 

Although  the  strong  and  electromagnetic  interactions  are  invariant  under  the 
C and  P operations,  it  came  as  quite  a surprise  to  physicists  to  find  that  the  weak 
interactions  were  not.  Order  was  restored  when  it  was  postulated  that  although 
the  weak  interactions  violated  the  C and  P symmetries  separately,  the  combined 
symmetry  CP  would  be  respected. 

1.3  Mixing  in  the  P°  — System 

Mixing  is  the  spontaneous  oscillation  of  a particle  into  its  anti-particle.  This 
behavior  can  happen  in  ground  state  electrically  neutral  mesons  with  net  fiavor,  such 
as  in  the  system  where  it  was  first  observed  [1],  for  example  through  “box” 

diagrams  like  the  one  shown  in  Fig.  1.1.  Other  systems  in  which  mixing  can  occur 
include  the  — P^,  and  P°  — P^  systems.  There  are  many  references 

describing  the  mixing  phenomenon,  particularly  in  the  and  P^  - P^  systems; 

some  good  examples  exist  [2,  3,  4,  5,  6].  In  this  section,  we  summarize  this  work  and 
see  what  it  can  tell  us  about  mixing  in  the  — P°  system. 

A similar  notation  can  be  used  to  describe  all  the  systems  in  which  mixing  can 
occur.  We  therefore  denote  the  eigenstates  of  flavor  as  P^{=K^,  D^,  B^)  and  P®. 

We  can  consider  a wave  function  made  of  an  arbitray  linear  combination  of  P°  and 
P^  at  time  t = 0,  and  express  that  wave  function  at  time  t = t as  a.  linear  combination 
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s,b  u,c,t  d 

Figure  1.1:  Box  Diagram  for  and  mixing. 


of  the  flavor  eigenstates  and  their  decay  products  |/) 


\m)  = + m\p°) + Y.  <^1^)  (1-5) 

/ 


where  the  |/)  are  the  states  into  which  and  can  decay.  As  the  system  evolves 
in  time,  and  P^  decay,  so  that  probability  migrates  into  the  decay  products  |/). 
We  write  the  full  Hamiltonian  H as  the  sum  of  two  pieces 

H = + (1.6) 


The  P[q  contains  the  strong  and  electromagnetic  interactions  and  has  |P^)  and  |F°) 
as  eigenstates.  The  contains  the  weak  interactions  responsible  for  the  decay  of 
P°)  and  IP*’)  states  and  their  transformations  among  each  other.  The  is  a small 
perturbation  on  Hq. 

Projecting  into  a space  spanned  by  |P°)  and  |P0),  the  evolution  of  |t/>(t))  is 
described  by  the  time-dependent  Schrodinger  equation 


. d 


a[t) 

j 


H 


eff 


Pii  Hu 

P21  H 


22 


' A(t)^ 


(1.7) 
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Assuming  CPT  invariance  gives 


where  Heir  is  the  projection  of  the  full  Hamiltonian  H into  the  subspace  spanned  by 
P®)  and  |P°).  We  have  denoted  the  diagonal  elements  of  M and  F as 


Fii  — T22  = Fo 


(1.10) 


which  stems  from  the  CPT  invariance.  Note  that  the  effective  Hamiltonian  is  not 
hermitian,  as  probability  can  leak  out  of  the  space  to  the  decay  products  |/). 

The  M is  called  the  mass  matrix,  and  F is  called  the  decay  matrix.  One  can 
show  by  performing  the  perturbation  theory  analysis  to  second  order  [6]  that  the 
elements  of  M and  F are  as  follows: 


M 


1-] 


= mpiy  + m.\]) + p E 

rup  — hf 


Tij  = 27t  ^ - Ef) 

I 


(1.11) 


where  V denotes  the  principal  value  prescription  is  applied  when  summing  over  the 
states  /,  avoiding  states  where  Ef  — mp.  In  addition,  we  used 


iFo|P°)  = mp|P") 


0' 


Po|PO)  = mp|po) 

Ho\f)  = Ef\f) 


(1.12) 
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The  values  of  the  matrix  elements  in  Equation  1.11  provide  insight  into  the 
physical  processes  described  in  mixing.  The  diagonal  terms  of  M are  dominated  by  the 
eigenvalue  mp  of  Ho,  which  makes  sense  as  the  strong  interactions  are  responsible  for 
binding  the  constituent  quarks  into  mesons.  The  off-diagonal  elements  of  M describe 
transitions  between  the  |P*^)  and  |P°)  (mixing!)  through  intermediate  states  |/). 
These  states  must  be  virtual  (off-shell),  as  the  principal  value  prescription  excludes 
the  on-shell  component  where  mp  = Ef.  The  term  (i\Hyj\j)  corresponds  to  a AS  = 2 
interaction.  Turning  to  the  decay  matrix  T,  we  can  see  that  diagonal  elements  arise 
from  the  allowed  decays  of  ^ / and  P°  with  the  states  / being  on-shell  due 
to  the  6{mp  — Ef).  Finally,  the  off-diagonal  elements  of  T describe  transitions  between 


P")  and  |P“)  (mixing!)  through  real,  on-shell  intermediate  states  /.  However,  / must 
be  a state  into  which  both  mesons  can  decay,  thus  describing  processes  P°  — > / ^ P^ 


and  P°  — ^ P^. 


We  can  solve  for  the  eigenstates  of  Heff . We  can  rewrite  Heff  as  a diagonal  piece 
and  an  off-diagonal  piece  (commonly  referred  to  as  the  mixing  matrix,  for  reasons  that 
soon  become  clear) 


We  note  that  eigenstates  of  the  off-diagonal  piece  are  also  eigenstates  of  Heff.  These 

eigenstates  are  the  mass  (also  known  as  physical)  eignestates;  these  are  the  states  we 
can  observe. 
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The  physical  eigenstates  are  written  as 


^ ^"’“0 

= 7M  ■ 

Solving  for  the  eigenvalues  of  Heff,  we  find 

(qV  _H2i  _ (M*,-iri,\ 

\pj  Hn  {Mu-iTu) 


(1.15) 


(1.16) 


and  the  eigenvalues  yUi,2  are 


Pi, 2 — Hu  ± (Hi2H2i)^  — Mo  ± (Hi2H2i)‘^  — M\^2  ~ (1-17) 

where 


= Mo±5R(Mi2//2i)^ 

= M„  ± Si 


and 


Ti,2  = roT2$5(Mi2M2i)^ 


= To  T 2^5 


(1.18) 


Note  that  the  signs  of  the  eigenvalues  do  not  have  to  correspond  to  the  signs 
denoting  the  eigenstates,  due  to  the  square  root.  The  mass  and  width  splittings 
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between  the  physical  eigenstates  are  then 


AM  = M_  - M+  = -25R  {HuH2iY 

Ar  = r_-r+  = 4S 


(1.19) 


For  completeness,  the  lifetimes  are  then 


n,2 


1 


ri,2 

1 


1 


ro  \ 1 T 


Ar 

2ro 


(1.20) 


and  the  physical  states  evolve  in  time  according  to 


^’i,2(«))  = = 0)). 


(1.21) 


Let’s  take  a moment  to  summarize  this  remarkable  result.  Particles  are  typically 
produced  as  eigenstates  of  flavor,  but  if  they  can  mix,  then  the  physical  eigenstates 
they  decay  as  are  not  the  same  as  the  states  in  which  they  were  produced!  The  physi- 
cal eigenstates  have  different  masses,  decay  widths,  and  lifetimes  than  the  eigenstates 
of  flavor.  Furthermore,  there  are  two  sources  of  mixing  (the  off-diagonal  terms  in  the 
mass  and  decay  matrices)  and  they  describe  very  different  physical  processes;  mixing 
through  either  real  or  virtual  intermediate  states,  respectively.  ^ 

The  limit  of  CP  invariance  for  H (and  therefore  Hen)  deserves  further  dis- 
cussion. In  this  limit,  M12  and  F12  are  real,  resulting  in  \q/p\  — 1.  The  physical 


®We  do  not  address  the  possibility  of  mixing  through  AP{=  S,  C,B)  =2  “superweak”  interac- 
tions. 
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eigenstates  |Pi,2)  become  the  eigenstates  of  CP,  |P±) 


p±)  = ^ ± \p“] 


(1.22) 


with 


CP|P_) 


(1.23) 


and  from  1.19  the  mass  and  width  splittings  between  the  physical  eigenstates  simplify 
to  \AM\  = 2Mi2  and  |Ar|  = 2rj2. 


1.4  Using  CP  Eigenstate  Decays  to  Probe  Mixing 

In  a P°  - P^  system  exhibiting  mixing,  the  relative  amount  of  constituent  flavor 
(or  equivalently,  CP)  eigenstates  in  an  initial  state  can  be  probed  with  decays  of  the 
physical  mass  eigenstates  into  CP  eigenstate  flnal  states. 

It  is  convenient  to  introduce  dimensionless  mixing  parameters  x and  y,  which 
are  deflned  as 


X = 


AM 

Ar 


y = 


2F 


0 


(1.24) 


We  have  already  seen  in  Equation  1.21  that  mixing  manifests  itself  partly  as  a 
difference  in  lifetime  between  the  flavor  eigenstates  \P^)  and  jP*^)  and  the  physical 
mass  eigenstates  |P±).  So,  how  to  observe  mixing?  We  can  interrogate  the  distri- 
bution of  times  in  the  decays  of  the  states  |P°(t)),  which  is  initially  pure  |P^),  and 


|P°(t)),  which  starts  at  t = 0 as  pure  |po).  Using  equations  1.15  and  1.21  and 
normalizing  p and  q such  that  |pp  + |g'p  = 1,  we  have 
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l^’V)>  = /+(i)l^’")  + -/-(^)|P"> 

P 

|p5(()>  = ^/_(i)|P»)  + U(t)\n  (1.25) 


where 


We  define  the  amplitudes  for  the  decays  of  the  fiavor  eigenstates  into  a final  state  / 
(P°  — )■  / and  — )■  /)  and  its  CP  conjugate  state  / (P°  — )■  / and  P'^  — >■  /)  as 


= (/|P|P°)  Af  = if\H\P0) 

Aj  = {7\H\P^)  Aj  = (7|Pp)  (1.27) 

then  the  decay  amplitudes  for  the  our  physical  eigenstates  are 

(/i/f|p»(())  = Mt)A,  + p-m, 

= ^Ut)A,  + U{t)A, 

Q. 

(f\H\P'‘{t))=U(t)Aj+^-f_{t)Aj 

J p J 

(7|/fF(()>  = ^f-it)Aj+U(t)Aj.  (1.28) 

For  convenience,  we  define  the  ratios 

pAf  ^ qAj 


(1.29) 
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then  the  time-dependent  decay  rates  are 


0 


oc 


l^/Pe 


2^— Tot 


(1  - |A/|2)  cos(AMt) 
-I-  (l  -H  |A/p)  cosh  ' 


- (A/)  sin  (AMt) 

+ 23fJ  (A/)  sinh 


R{P%t)^f)<x\{f\H\P0{t))\^ 


OC 


p 

2 r 

Q 

(|A/p  — l)  cos  (AMt) 
-I-  (l  -I-  |A/p)  cosh  ' 


+ 2^  (A/)  sin  (AMt) 

Art 


-I-  25?  (A/)  sinh 


and  for  decays  into  / final  states 


i?(P«(t)^7)cx|(/|//|P°(t))| 


0 


OC 


1^ 


—\2p—rot 

f\  e 


Q 

P 


(|Ajp  — l)  cos  (AMt) 

Art 


-I-  (l  -H  |Ajp)  cosh  r 

-h  2^  (Aj)  sin  (AMt) 

, Art 

H-  23?  (Aj)  sinh 


(1.30) 
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R(P“{t)  ^ /)  a K/|J/|P»(())|^ 


OC 


(l  — |Ajp)  cos  (AMt) 

/i  I X i2\  1 ' 

+ (l  + |Ajf ) cosh 

— 2S5  (Aj)  sin  (AMt) 
+ 2U  (Xj)  sinh  ' 


(1.31) 


It  is  obvious  from  the  the  equations  above  that  the  distribution  of  decay  times  gener- 
ally deviates  from  a pure  exponential,  and  the  deviation  is  expressed  in  terms  of  the 
quantities  which  characterize  mixing,  AM  and  AF.  If  AM  and  AF  are  both  zero  (no 
mixing),  the  distribution  of  decay  times  returns  to  a simple  exponential. 


If  we  return  to  the  limit  of  CP  invariance,  then 
tions  1.30  (and  1.31)  can  be  combined,  resulting  in 


E 

Q 


= 1.  The  rates  in  equa- 


R{P\t)  ^f)  + R{P0  ^ /)  OC 


RiP^{t)  ^ /)  + i?(P0  ^ /)  OC  \Aj\h-^°^ 


(l  + |A/|^)cosh 
+ 2U{Xf)smh(^^\ 
(1  + \Xj\^)  cosh 

/ AFt 

+ 23?  (Aj)  sinh 


(1.32) 


One  can  see  that  by  studying  the  combined  time-dependent  rates,  we  are  sensitive 
only  to  AF.  Similarly,  by  studying  the  difference  in  rates,  we  are  sensitive  to  AM. 
Further,  if  / is  a CP  eigenstate  (|/)  = CP\f)  = ricp{f)\f)),  then  the  distribution  of 
times  (in  the  “combined”  case)  can  be  expressed  as 
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R{P^{t)  ^ /)  + R{P^  ^ /)  oc  2\Af\‘^e-^°^ 


cosh 


+ ?7cp(/)sinh 


(1.33) 


where  the  subscript  CP  has  been  added  to  y to  denote  that  CP  invariance  is  assumed. 
We  note  this  rate  is  again  a simple  exponential! 

Also  in  the  CP  conserving  limit,  we  can  regard  decays  to  non-CP  eigenstates 
as  equal  admixtures  of  CP  even  and  odd  components;  call  them  UP.  We  can  then 
rewrite  ycp  as 


ycp  = - 1 CP  Even 

TCP+ 

ycp  = l-  CP  Odd 

TCP- 


(1.34) 


It  is  important  to  note  that  in  the  CP  conserving  limit, 


P+)(|P_))  corresponds  to 


So,  in  the  limit  of  CP  invariance,  mixing  due  to  the  difference  in  the  widths 
of  the  physical  eigenstates  can  be  seen  as  the  difference  in  lifetimes  between  a CP 
eigenstate  final  state  and  another  final  state  which  is  not  a CP  eigenstate  (or  has 
different  eigenvalue  t]cp)-  This  is  our  analysis.  We  compare  the  lifetimes  of  CP  even 
states  KK  and  tttt  to  that  of  CP  state  Ktt. 

1.5  Theoretical  Expectations  for  ycp  and 
Current  Experimental  Limits 

A large  collection  of  theoretical  mixing  predictions  are  summarized  in 

reference  [7],  covering  a time  range  of  nearly  20  years.  There  are  several  things  to 
note  about  this  these  predictions.  First,  mixing  predictions  for  x and  y tend  to  be 


18 


small,  and  cover  many  orders  of  magnitude  10“^  - for  either  quantity).  We 
expect  both  of  these  parameters  to  be  small  in  the  Standard  Model;  x should  be  small 
because  all  the  quarks  in  the  box  diagram  have  masses  that  are  basically  degenerate 
on  the  scale  of  Mw,  and  y should  be  small  because  long  distance  effects  are  expected 
to  be  negligible  [8].  Second,  the  upper  limit  of  all  mixing  predictions  is  within  reach 
of  current  experiments.  Third,  non-Standard  Model  physics  tends  to  enhance  the 
predictions  for  x.  This  is  because  such  models  tend  to  allow  heavy  particles  into  the 
box  diagram.  The  prospects  for  y,  on  the  other  hand,  do  not  improve  with  such  the 
additions  to  the  particle  zoo  [9]. 

A plot  of  current  experimental  limits  on  y in  the  CP  conserving  limit  (ycp) 
from  studies  of  D meson  decays  into  CP  eigestates  is  shown  in  Figure  1.2.  These 
results  are  documented  in  references  [8,  10,  11].  The  BELLE,  FOCUS,  the  E791  KK 
analyses  are  of  the  same  style  as  the  one  presented  here  (though  we  use  the  tttt  decay 
channel,  as  well).  Therefore,  these  results  are  directly  comparable  to  the  one  we 
present  later.  The  tightest  limits  shown  on  Figure  1.2  come  from  FOCUS.  Expressed 
as  a 95%  confidence  level  interval,  FOCUS  finds: 

0.0034  < ycp  < 0.0650 

This  is  not  quite  consistent  with  zero,  which  is  very  interesting.  The  FOCUS  result 
is  the  first  result  to  be  even  mildly  inconsistent  with  zero. 

Other  analyses  have  also  provided  insight  on  mixing.  The  CLEO 

“wrong  sign  D”  analysis  [12]  is  not  sensitive  to  a:  or  y by  themselves,  but  to  the  total 
mixing  rate: 


Rmix  = y/x'^  -\-y‘^  = y/ X'^  + ?/'2 
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where  x'  and  y'  are  just  x and  y rotated  by  <5: 

x'  = xcos5  + ysinS 
y'  = ycos5  — xsin^ 

and  is  a stong  phase  diference  between  the  right-  and  wrong-sign  A'tt  decay  ampli- 
tudes. At  the  95%  confidence  level,  this  result  is; 

(1/2)  x'^  < 0.00041 
-0.058  <y'  < 0.010 

with  no  assumption  about  CP  violation. 

The  E791  Klv  analysis  [13]  is  sensitive  only  to  Rmix\  the  limits  derived  from 
this  analysis  carves  out  a semi-circle  in  the  xy  plane.  They  find 


Rmix  < 0.0050 


at  the  90%  confidence  level. 

To  summarize,  with  current  experimental  sensitivities  (particularly  from  CLEO 
and  FOCUS)  in  the  few  percent  range,  it  may  be  that  a mixing  signal  is  just  around 
the  corner.  If  a signal  is  seen  above  the  10~^  level,  it  would  favor  a “new  physics” 
explanation. 
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CP  Eigenstate  Summary 


ycp(%) 


Figure  1.2;  Experimental  limits  on  ycp  from  CP  eigenstate  decays  of  D 

mesons.  The  error  bars  shown  are  Icr  errors  (the  statistical  and 
systematic  errors  added  in  quadrature). 


CHAPTER  2 

EXPERIMENTAL  APPARATUS 


The  D mesons  we  want  to  study  are  produced  only  in  a high  energy  environment 
and  decay  very  quickly  into  daughter  particles  which  fly  away  at  relativistic  speeds.  In 
order  to  study  such  particles,  an  environment  capable  of  producing  and  detecting  such 
particles  and  detecting  must  be  created.  There  are  tremendous  technical  challenges 
in  creating  such  an  apparatus  and  in  analyzing  the  data  which  comes  out  of  it.  The 
solutions  to  these  challenges  are  many,  and  as  one  might  guess  they  have  evolved  over 
time  as  technology  has  progressed.  As  documenting  the  complete  history  of  particle 
physics  techniques  is  clearly  beyond  the  scope  of  this  work,  we  instead  focus  in  this 
chapter  on  the  apparatus  used  for  this  analysis,  namely  the  CESR  e‘^e~  collider  and 
the  CLEO  detector. 

2.1  The  Cornell  Electron  Storage  Ring 

The  Cornell  Electron  Storage  Ring  (CESR)  is  a circular  electron-positron  col- 
lider located  at  Cornell  University’s  Wilson  Laboratory  in  Ithaca,  New  York.  The 
layout  of  Wilson  Laboratory’s  accelerator  facility  is  shown  in  Fig.  2.1. 

CESR  is  staffed  by  accelerator  physicists  and  students  primarily  involved  in 
the  development,  operation,  maintenance,  and  construction  of  accelerator  machinery. 
CESR  itself  has  been  around  for  quite  a long  time;  it  has  provided  e'''e“  collisions  and 
synchrotron  radiation  for  several  experiments  since  it  came  online  in  1979.  During  the 
decade  of  the  1990s,  which  is  the  time  scale  relevant  for  this  analysis,  e'^e“  collisions 
occurred  only  in  the  south  interaction  region  (at  the  bottom  of  Fig.  2.1),  where  the 
CLEO  II  detector  is  depicted. 
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O Positron  Bunch  - Clockwise 
# Electron  Bunch  - Counter  Clockwise 


Figure  2.1:  Schematic  of  the  Wilson  Laboratory  accelerator  facility.  CESR 

and  the  synchrotron  are  in  the  same  tunnel,  about  40  feet  beneath 
Cornell  University’s  Alumni  Fields.  The  circumference  of  CESR 
is  768  meters. 
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Running  CESR  in  “HEP”  (high  energy  physics)  mode  involves  creating,  accel- 
erating, storing,  and  colliding  beams  of  electrons  and  positrons  so  that  the  resident 
experiments  can  take  data.  These  processes  are  accomplished  serially  and  then  re- 
peated for  several  days  or  even  weeks  at  a time.  “HEP”  mode  is  normally  interrupted 
by  scheduled  periods  of  “machine  studies”  for  study  and  experimentation  with  new 
running  conditions  (there  is  always  something  to  learn  and  improve)  as  well  as  regu- 
lar maintenance  of  the  accelerator  components;  such  down  times  typically  last  a few 
days.  More  extensive  maintenance  and  upgrades  to  the  accelerator  requiring  down 
times  lasting  a few  to  many  weeks  are  less  common,  happening  perhaps  a couple 
times  per  year. 

The  act  of  creating,  accelerating  to  operating  energy,  and  storing  of  electron  and 
positron  beams  is  known  as  a “fill”.  Normally,  the  fill  process  entails  the  replenishing 
(or  creation)  of  electron  and  positron  currents  lost  to  some  nominal  amount.  Once 
the  beams  are  then  prepared  for  collision,  or  “tuned”,  collisions  can  begin.  When  in 
HEP  mode,  CESR  provides  colliding  beams  for  a predetermined  length  of  time;  this 
is  what  is  known  as  a “run”.  Over  the  years,  the  length  of  a run  has  varied  somewhat, 
but  has  generally  been  on  order  of  one  hour.  When  a run  is  ended,  the  process  begins 
again.  During  the  1990s,  there  were  over  31,000  HEP  runs. 

As  Fig.  2.1  suggests,  CESR  is  composed  of  several  major  components,  each  of 
which  performs  a specific  function  in  the  arduous  task  of  creating  collisions. 
These  components  and  their  purpose  are  described  in  the  subsections  below. 

2.1.1  Linear  Accelerator 

In  order  for  e''‘e“  collisions  to  occur,  electrons  and  positrons  must  somehow  be 
produced  in  the  first  place.  This  is  the  job  of  a special  30  meter  section  of  vacuum 
pipe  called  the  Linear  Accelerator  (LINAC),  and  the  process  is  called  “filling”  the 
synchrotron.  The  fill  procedure  is  sequential:  positrons  first,  then  electrons. 
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The  positron  filling  procedure  is  as  follows.  First,  electrons  are  created  by 
boiling  them  off  of  a filament  wire  located  at  the  back  of  the  LINAC  (the  electron 
gun),  just  like  inside  the  picture  tube  of  a television.  These  electrons  are  accelerated 
by  electromagnetic  fields  in  a series  of  Radio  Frequency  Acceleration  Cavities  (RF 
Cavities)  towards  a tungsten  target  located  about  halfway  down  the  LINAC.  The 
collision  between  the  140  MeV  electrons  and  the  target  creates  a spray  of  electrons, 
positrons,  and  photons.  The  electrons  are  cleared  away  with  magnetic  fields  and  the 
positrons  are  accelerated  down  the  remainder  of  the  LINAC,  ultimately  reaching  an 
energy  of  about  200  MeV . At  the  end  of  the  LINAC,  the  positrons  are  introduced  into 
the  synchrotron.  The  electron  filling  procedure  is  identical,  except  that  the  tungsten 
target  is  moved  out  of  the  way;  the  conversion  process  is  not  necessary,  as  we  already 
have  electrons  from  the  electron  gun.  The  total  filling  time  for  electrons  and  positrons 
in  normal  running  conditions  is  about  10  minutes. 

We  skipped  an  important  detail.  The  astute  reader  probably  noticed  the  tremen- 
dous bulbs  on  the  CESR  ring  in  Fig.  2.1  labelled  “Positron  Bunch”  and  “Electron 
Bunch”.  No,  they  are  not  really  that  big;  the  bunches  of  positrons  and  electrons  fit 
quite  comfortably  inside  of  all  of  CESR’s  vacuum  pipes.  The  electron  and  positron 
beams  are  not  continuous;  instead  they  are  grouped  in  bunches.  The  creation  of  these 
bunches  begins  in  the  LINAC’s  “prebuncher”  component.  The  prebuncher  is  a short 
series  of  cavities  positioned  immediately  outside  the  electron  gun.  When  filling  the 
electrons  and  positrons,  the  electron  gun  is  fired  repeatedly;  each  pulse  resulting  in  a 
stream  of  electrons  approximately  3 feet  long.  The  electron  pulses  are  compressed  in 
the  prebuncher  with  electromagnetic  fields  to  the  width  of  a dime. 

2.1.2  Synchrotron 

The  synchrotron  is  a circular  accelerator  in  which  electrons  and  positrons  from 
the  LINAC  are  accelerated  to  operating  energy.  The  electrons  and  positrons  and 
travel  inside  a vacuum  pipe  in  opposite  directions  (viewed  from  above,  electrons  go 
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counterclockwise,  and  the  positrons  clockwise)  in  circular  orbits.  Steering  is  accom- 
plished via  magnetic  field,  and  acceleration  by  radio  frequency  electromagnetic  field. 
As  the  particles’  speed  increases,  these  fields  must  be  adjusted  in  a synchronized 
manner  to  keep  the  radius  of  the  orbit  constant. 

For  most  of  the  circumference  of  the  ring  the  particles  travel  in  vacuum  between 
the  poles  of  bending  magnets,  arranged  in  192  3-meter  long  sections.  At  four  places 
in  the  ring  the  beam  passes  through  a 3-meter  long  linear  accelerator.  After  approx- 
imately 4000  laps  around  the  synchrotron  (about  1/lOOth  of  a second),  electron  and 
positron  bunches  have  reached  operating  energy,  about  5.3  GeV  each,  and  are  ready 
for  injection  into  the  storage  ring. 

2.1.3  Storage  Ring 

The  positron  and  electron  beams  are  transferred  into  the  storage  ring  (CESR) 
via  transfer  lines  that  connect  the  synchrotron  to  CESR.  This  process  is  known  as 
“injection”.  There  is  one  transfer  line  for  positrons  (“west”)  and  one  for  electrons 
(“east”).  At  this  point,  the  electron  and  positron  buches  are  travelling  at  speeds  over 
99.9999995%  of  the  speed  of  light.  The  purpose  of  the  storage  ring  is  to  keep  the 
beams  of  electron  and  positron  bunches  constant  and  at  the  same  time  arrange  for  a 
high  rate  of  e'*"e“  collisions  in  the  interaction  region. 

The  storage  ring  uses  the  same  guide  field  principle  as  the  synchrotron;  that 
is,  particles  travel  in  a (roughly)  circular  orbit  in  vacuum  under  the  influence  of  a 
magnetic  field.  Each  electron  or  positron  spends  less  than  one-hundredth  of  a second 
in  the  synchrotron,  but  once  it  is  transferred  to  the  storage  ring  it  must  coast  there 
for  up  to  several  hours.  This  puts  very  strict  demands  on  the  precision  of  the  ring  of 
magnets  and  on  the  quality  of  the  vacuum.  The  electrons  and  positrons  can  strike 
and  scatter  off  of  any  gas  molecules  they  encounter  in  the  beampipe.  If  there  are 
too  many  gas  molecules,  this  ultimately  leads  to  loss  of  the  beams.  Therefore,  the 
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residual  air  in  the  path  of  the  beam  particles  must  be  kept  about  one  hundred  million 
times  less  dense  than  in  the  normal  atmosphere. 

The  orbits  of  the  electrons  and  positrons  bring  them  into  collision  only  at  the 
interaction  region  (at  the  center  of  the  CLEO  detector).  As  mentioned  above,  a 
magnetic  fields  are  used  to  make  the  particle  orbits  roughly  circular;  a total  of  86 
dipole-type  magnets  are  used  to  accomplish  this.  Another  106  quadrupole  magnets 
are  used  to  keep  the  beams  focused,  as  we  wish  the  particles  to  remain  in  tight 
bunches.  A small  number  of  sextupole  and  octupole  magnets  make  subtle  corrections 
to  the  orbits.  The  sequence  of  dipoles  and  quadrupoles  repeats  itself  around  the  768 
meter  circumference  of  the  ring. 

The  particles  coasting  around  the  storage  ring  radiate  a beam  of  X-rays  like  a 
headlight  on  a train,  losing  about  1.2  MeV  per  turn  around  the  ring.  The  power 
level  of  this  “synchrotron  radiation”  is  about  half  a megawatt.  It  strikes  the  vacuum 
chamber  in  a narrow  stripe,  depositing  energy  on  the  vacuum  chamber  wall.  The  heat 
generated  must  be  carried  away  by  water  circulating  through  a channel  in  the  vacuum 
chamber  wall.  The  lost  energy  is  restored  to  the  circulating  beam  by  radio-frequency 
(RF)  cavities  similar  to  those  in  the  synchrotron  and  linac,  although  the  beam  is  not 
actually  accelerated.  The  RF  system  operates  at  a frequency  of  500  million  cycles  per 
second  (MHz)  and  is  powered  by  a system  rather  like  a television  transmitter.  The 
intense  synchrotron  radiation  is  not  simply  a useless  byproduct  of  acceleration;  it  is 
used  for  experiments  in  the  CHESS  area  (to  the  left  and  right  of  CLEO  in  Fig  2.1). 

We  have  previously  stated  that  electrons  and  positrons  travel  around  the  storage 
ring  in  bunches,  as  opposed  to  a continuous  stream.  In  fact,  there  is  actually  finer 
structure  than  this.  The  bunch  structure  is  more  appropriately  described  as  equally 
spaced  trains  with  up  to  five  cars  each  that  can  be  populated  with  electrons  and 
positrons.  In  each  train  (or  “bunch-train”),  the  cars  are  closely  spaced,  about  14  ns 
apart.  This  fact  requires  that  when  the  beams  are  ultimately  crossed  in  the  interaction 
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region,  that  they  cross  at  an  angle  in  order  to  prevent  “parasitic”  crossings  away  from 
the  nominal  interaction  point.  In  the  1990s,  CESR  ran  with  nine  trains  of  electrons 
and  positrons  and  up  to  three  cars  populated  with  approximately  around  1.5  x 10^^ 
particles  each  (mid-1990s  running).  At  5.3  GeV^  the  cars  are  circling  the  storage 
ring  at  396  kHz.  It  should  be  noted  that  CESR  developed  the  idea  of  bunch-train 
operations  and  was  the  first  accelerator  to  use  them  in  production.  It  is  one  of  many 

collision  rate  enhancement  techniques  CESR  has  come  up  with  which  has  been 
adopted  at  other  laboratories,  including  CERN,  SLAC,  and  KEK. 

Certainly,  we  do  not  want  the  beams  to  collide  anywhere  besides  the  interaction 
region  (where  the  detector  is),  so  the  orbits  are  not  perfect  circles.  Large  capacitors 
called  “horizontal  separators”  are  used  to  knock  the  bunches  slightly  out  of  the  circle, 
the  species  moving  in  opposite  directions.  In  the  end,  the  orbits  are  shaped  like  a 
pretzel. 

The  behavior  of  collections  of  charged  particles  in  a storage  ring  is  a complicated 
field  of  its  own,  whose  details  are  certainly  beyond  the  scope  of  this  work.  Suffice  it 
to  say  that  we  have  not  discussed  many  of  the  important  effects  which  complicate  the 
operation  of  a real-world  storage  ring,  or  the  technically  challenging  solutions  which 
must  be  implemented  to  make  it  all  work.  To  that  end,  the  CESR  staff  has  constructed 
a system  with  sophisticated  computer  control,  monitoring,  and  feedback  systems 
which  communicate  with  the  magnets  and  the  RE  system  in  all  of  the  accelerator 
components. 

2. 1.3.1  Interaction  Region 

The  interaction  region  (IR)  is  where  the  electron  and  positron  bunches  pass 
through  each  other,  which  they  do  about  3.6  million  times  per  second.  The  IR  is  in  a 
special  segment  of  the  accelerator  vacuum  system;  it  is  the  part  which  passes  through 
the  heart  of  the  CLEO  detector.  Indeed,  this  beam  pipe  is  as  much  a part  of  CLEO 
as  it  is  of  CESR.  This  is  because  the  particles  created  in  an  e'''e“  annihilation  must 
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Figure  2.2:  The  CLEO  II. V beam  pipe. 


Figure  2.3:  The  CLEO  II. V beam  pipe  surrounded  by  the  Silicon  Vertex 

Detector  and  Vertex  Detector. 

first  pass  through  the  walls  of  this  segment  of  the  beam  line  before  encountering  the 
instrumented  CLEO  detector  material.  The  pipe  is  therefore  made  of  a low  mass 
material  and  is  very  thin  in  order  to  reduce  scattering.  If  this  were  not  the  case, 
charged  particles  emerging  from  the  e+e“  interaction  point  and  subsequently  tracked 
in  the  CLEO  tracking  system  would  not  point  back  to  the  place  they  began.  There 
were  two  beam  pipes  used  during  the  CLEO  II’s  lifetime.  The  first  version  was  used 
from  1989  to  1995;  it  had  an  inner  radius  of  3.5  centimeters  and  was  0.5  millimeters 
thick.  The  second  was  installed  as  part  of  “Phase  11”  of  the  CESR/CLEO  Upgrade 
Project  [18]  in  the  summer  of  1995.  The  inner  radius  of  the  beam  pipe  was  reduced 
to  1.875  centimeters,  and  the  thickness  was  approximately  1 millimeter.  The  second 
beam  pipe  was  actually  double- walled,  so  it  isn’t  actually  twice  as  thick  as  the  previous 
version.  Between  the  walls,  water  was  used  for  cooling  in  order  to  compensate  for 
the  extra  heat  resulting  from  higher  beam  currents.  Both  beam  pipes  were  made  of 
beryllium.  Figures  2.2  and  2.3  show  pictures  of  the  second  version  of  the  beam  pipe. 
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If  the  accelerator  is  performing  well,  a high  collision  rate  is  the  result.  Indeed, 
the  ability  to  have  a high  collision  rate  is  crucial  for  the  success  of  an  accelerator  and 
the  experiments  it  serves.  The  figure  of  merit  is  then  the  number  of  possible  collisions 
per  second  per  unit  area;  this  is  called  the  luminosity.  The  higher  the  luminosity,  the 
larger  the  rate  for  collisions.  And  the  larger  the  number  of  collisions,  the  larger 
the  data  samples  we  have  to  study.  The  luminosity  is  given  as 
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where  / is  the  frequency  of  revolution  for  each  train,  n is  the  number  of  populated 
cars  in  each  train  for  each  particle  species,  A is  the  cross-sectional  area  of  the  cars, 
and  Ne+  and  N^-  are  the  numbers  of  positrons  and  electrons  per  car,  respectively. 

In  order  to  maximize  the  luminosity,  the  beams  are  focused  to  be  as  small  as 
possible  at  the  IR.  Typically,  each  car  is  about  2 cm  long  (along  the  beamline),  300  /xm 
wide  (in  the  plane  of  the  ring),  and  8 /xm  high  (normal  to  the  ring).  Over  time,  CESR 
has  found  ways  to  increase  the  luminosity  while  retaining  stable  running  conditions 
by  adding  more  current  (increase  N^+  and  N^-)  and  filling  more  cars  (n).  Fig.  2.4 
shows  the  CESR  luminosity  hisory,  as  seen  by  the  CLEO  II. V experiment. 

During  the  1990s,  CESR  held  and  constantly  improved  upon  the  world  record 
for  luminosity  until  the  February  1999  shutdown  for  the  next  incremental  upgrade 
to  the  machine  and  the  CLEO  detector.  The  highest  luminostiy  achieved  by  CESR 
before  the  shutdown  was  8.0  x 10^^ cm~^ sec~^ , exceeding  the  Phase  II  design  goal  by 
33%.  The  excellent  performance  of  CESR  during  this  time  allowed  the  accumulation 
of  far  more  data  than  at  any  other  experiment.  There  are  many  accelerators  of  higher 
energy  than  CESR,  but  there  are  none  which  even  came  close  to  the  luminosity  CESR 
delivered  in  this  time.  All  this  on  a tight  budget  to  boot.  In  the  future,  however, 
CESR  will  have  to  fight  hard  for  the  luminosity  crown  with  new  accelerators  online  at 
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Figure  2.4:  Daily  luminosity  delivered  to  the  CLEO  II. V experiment. 
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Figure  2.5:  Cross  section  versus  Mass. 


SLAC  and  KEK.  At  the  time  of  this  writing,  the  new  PEP-II  storage  ring  at  SLAC 
has  shattered  the  CESR  record  by  attaining  luminosity  of  1.43  x 10^^cm“^sec“F 
Meanwhile,  the  Phase  III  upgrade  of  CESR  has  been  installed  is  being  commissioned. 

2.2  The  Upsilon  System 

Why  is  the  CESR  beam  energy  5.3  GeVl  As  one  might  imagine,  there  is  a very 
good  reason  for  this. 

As  originally  designed,  CESR  could  run  at  center  of  mass  energies  as  high  as 
16  GeV.  However,  at  a fixed  target  experiment  at  Fermilab  in  1977  the  T family 
of  resonances  was  discovered.  These  resonances  are  shown  in  Fig.  2.5.  Though  the 
Fermilab  experiment  was  unable  to  resolve  the  members  of  this  new  family  as  separate 
entities,  they  were  immediately  believed  to  be  the  bound  states  of  a new  flavor, 
bottom,  just  as  the  observation  of  the  J/'ip  signalled  the  discovery  of  charm.  CESR’s 
first  achievement  was  the  resolution  of  the  first  three  states  in  T family;  the  T(1S), 
T(2S),  and  T(3S).  Discovery  of  the  T(4S)  followed  shortly  in  1980.  The  T(4S)  turns 
out  to  be  a very  interesting  state  indeed. 


32 


In  Fig.  2.5,  one  can  tell  that  there  is  an  obvious  difference  in  the  widths  between 
the  T(4S)  and  the  T(1S),  T(2S),  or  T(3S).  This  implies  that  the  T(4S)  has  decay 
channels  open  to  it  that  the  others  do  not  have.  These  decay  modes  are  and 

B^B^,  which  the  T(4S)  is  just  massive  enough  to  produce.  T(4S)  decays  nearly  100% 
of  the  time  to  B~^B~  and  B^B^  in  approximately  equal  proportions.  This  makes  the 
T(4S)  an  ideal  laboratory  to  study  the  characteristics  of  the  B meson,  and  is  why 
CESR  runs  with  the  beam  energies  it  does. 

Fig.  2.5  is  interesting  for  another  reason.  One  can  see  that  the  T resonances 
are  all  sitting  on  a considerable  amount  of  hadronic  “background” . This  background 
is  composed  primarily  of  e'*'e“  — >•  99,  where  g is  a light  quark  {udsc).  We  call  this 
process  “continuum”.  There  is  also  some  contribution  from  e+e“  r+r“,  where 
one  of  the  r leptons  decays  into  hadronic  products.  Hadronic  data  collected  at  the 
T(4S)  os  dominated  by  continuum,  where  the  cross  section  outnumbers  that  of  BB 
production  by  about  4:1. 

CESR  spends  about  2/3  of  its  time  running  on  the  T(4S)  resonance,  with  the 
remainder  of  the  time  spent  just  below  it.  This  is  done  in  order  to  collect  large 
amounts  of  data  we  know  contains  no  B mesons.  The  study  of  the  continuum  is  very 
important  to  B meson  analyses.  A B meson  physicist  ultimately  wishes  to  remove 
the  continuum  as  much  as  possible;  the  continuum-only  data  is  an  independent  data 
sample  which  allows  them  to  study  how  to  best  go  about  that. 

But  in  this  analysis,  the  tables  are  reversed.  We  want  to  measure  the  lifetime. 
To  do  so,  we  need  to  know  where  it  starts  with  high  precision.  For  D mesons  which 
are  decay  products  of  B mesons,  we  cannot  measure  the  D production  point  {B  decay 
vertex)  without  reconstructing  the  B meson,  which  simply  cannot  be  done  with  high 
efficiency.  Therefore  it  is  the  continuum  which  is  interesting  to  us,  and  the  BB  events 
are  a background  we  must  find  a way  to  remove  from  our  data  sample. 
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2.3  The  CLEO  Experiment 

In  an  e+e"  annihilation  at  CESR  energies,  the  vast  majority  of  the  time  the  e+e” 
interact  electromagnetically,  creating  a photon  This  interaction,  like  all  others, 
must  conserve  energy  and  momentum.  The  produced  photon  then  has  energy  which 
is  the  sum  of  the  energies  of  the  e'*'  and  e~  (in  the  neighborhood  of  10.58  GeV),  and 
a very  tiny  momentum,  as  the  incident  and  e~  momenta  nearly  cancel  each  other 
out.  The  energy  of  a particle  is  related  to  its  rest  mass  and  momentum  through  the 
relativistic  energy-momentum  equation,  written  as 

E = + m?  (2.2) 

This  means  that  the  photon  produced  from  an  e'^'e"  collision  at  CESR  has  mass.  This 
is  a common  situation  in  an  e“''e~  collider;  indeed  this  fact  is  exploited  in  the  design 
of  such  facilities.  Such  a photon,  written  as  7*,  is  referred  to  as  “off  mass  shell”  or  as 
a “virtual”  photon  because  real  photons  have  mass  of  zero  (see  Table  1.3).  Virtual 
photons  are  unstable,  and  decay  immediately  into  “real”  (or  “on  shell”)  daughters. 
It  is  at  this  point  we  enter  the  world  of  particle  detectors. 

Particle  detectors  everywhere  have  at  least  one  thing  in  common;  they  are  made 
of  matter.  As  we  described  in  Chapter  1,  matter  can  interact  with  other  matter. 
Particle  detectors  are  sensitive  to  such  interactions  and  are  instrumented  to  record 
information  about  the  interaction.  Using  this  information,  humans  can  determine 
some  of  the  properties  of  the  incident  particle,  such  as  its  energy,  momentum,  and 
mass. 


^Recall  from  Table  1.2  that  e+  and  e~  both  have  weak  isospin,  and  can  interact  weakly.  It  is 
certainly  true  that,  even  though  CESR’s  energy  scale  is  small  compared  to  the  mass  of  the  an 
e+e“  annihilation  could  result  in  an  off-shell  Z°.  However,  the  ratio  of  the  electromagnetic  to  weak 
coupling  constants  is  large,  so  electromagnetic  interactions  dominate 
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There  are  many  ways  in  which  the  7*  can  decay.  In  this  analysis,  the  process 
of  interest  is  e''‘e“  ^ 7*  ^ hadrons  where  one  or  more  of  the  hadrons  or  its  decay 
products  is  a (or  D^)  meson.  The  D mesons  themselves  are  unstable,  and  quickly 
decay  into  lighter  particles.  It  is  these  particles  that  we  detect. 

In  this  analysis,  we  use  the  CLEO  II. V detector  to  find  the  daughters  of  the  D 
mesons.  In  the  following  sections,  we  discuss  some  of  the  particle  detection  schemes 
and  techniques  implemented  in  the  CLEO  detectors,  and  how  the  raw  detector  data 
is  transformed  into  measurements  of  particle  energy,  momenta,  and  trajectories. 

2.3.1  The  CLEO  II  Detector 

The  CLEO  II  detector  is  a versatile,  general  purpose  detector  of  cylindrical 
geometry  with  excellent  charged  particle  and  photon  detection  capabilities.  It  is 
thoroughly  described  in  reference  [14].  As  implied  by  the  name,  CLEO  II  is  not  the 
only  CLEO  detector  which  has  existed.  CLEO  II  was  preceded  by  CLEO  I.V  and  the 
original  CLEO  detectors.  Compared  with  the  previous  versions  of  CLEO,  CLEO  II 
was  a major  upgrade,  with  many  detector  elements  replace  with  new  ones. 

A side  view  of  CLEO  II  is  shown  in  Fig.  2.6,  and  a close  up  view  of  one  quadrant 
is  shown  in  Fig.  2.7.  As  one  can  see,  CLEO  II  is  a composite  of  many  detector 
elements.  These  subdetectors  are  typically  arranged  as  concentric  cylinders  centered 
around  the  beam  pipe.  The  entire  detector  is  approximately  cube  shaped,  with  one 
side  measuring  about  6 meters,  and  weights  over  1000  tons.  CLEO  II  operated  in 
this  configuration  from  1989  to  1995. 


^The  hadronization  mechanisms  describing  the  production  of  hadrons  from  the  7*  will  not  be 
discussed  in  this  document.  Suffice  to  say  that  e"'‘e“  — t 7*  — > hadrons  is  just  one  of  the  many 
processes  possible  in  e"*'e“  annihilation,  and  is  the  process  of  interest  in  this  analysis. 
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Figure  2.6:  A side  view  of  the  CLEO  II  detector. 
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Figure  2.7:  A side  view  of  a quadrant  of  the  CLEO  II  detector. 


It  will  be  useful  at  this  point  to  define  the  coordinate  system  relavent  to  the 
CLEO  detector.  The  x direction  points  south  of  the  ring  at  the  interaction  point,  y 
points  upwards,  and  2 is  the  direction  of  the  positrons  (west).  The  azimuthal  angle  (f) 
is  defined  in  the  canonical  way,  with  (f)  = 0 along  the  +x  direction.  The  polar  angle 
9 is  measured  from  the  positive  z axis. 

2. 3. 1.1  Superconducting  Coil 

All  the  CLEO  detector  subsystems  except  for  the  muon  chambers  are  located 
inside  a superconducting  coil.  The  purpose  of  the  coil  is  to  provide  the  1.5  Tesla 
magnetic  field  used  to  bend  the  paths  of  charged  particles  in  the  tracking  system. 
The  magnetic  field  points  in  the  -2  direction  (east).  By  measuring  how  much  a 
charged  particle  bends,  experimenters  can  measure  the  momentum  of  the  particle. 
The  helium  reservoir  shown  in  Fig.  2.6  at  the  top  of  CLEO  II  contains  the  liquid 
helium  required  to  keep  the  coil  in  a superconducting  state. 
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2. 3. 1.2  Charged  Particle  Tracking  System 

After  particles  from  the  interaction  point  pass  through  the  beam  pipe,  they 
begin  to  encounter  the  active  detector  elements  of  the  tracking  system.  In  CLEO  II, 
there  are  three  main  components  of  the  tracking  system,  all  designed  to  measure  a 
charged  particle’s  position  at  various  points  along  its  trajectory;  the  tracking  system 
is  not  sensitive  to  particles  with  no  electric  charge.  The  tracking  system’s  components 
all  fall  into  a category  of  tracking  devices  known  as  drift  chambers.  As  a result,  they 
all  have  the  same  principle  behind  their  operation. 

A drift  chamber  is  a volume  of  gas  (typically  a mixture  of  organic  and  inert 
gases)  with  an  array  of  anode  wires  and  cathode  wires^  or  metallic  surfaces  running 
through  it.  The  anodes  are  kept  at  high  (positive)  voltage,  which  provides  an  electric 
field  throughout  the  chamber.  The  cathodes  are  held  at  ground,  and  shape  the  electric 
field  so  that  fields  from  neighboring  anode  wires  do  not  interfere  with  each  other.  The 
wires  are  often  arranged  in  “drift  cells”  centered  on  the  anode.  Typically,  there  are 
several  to  many  layers  of  such  cells  in  a drift  chamber,  with  each  layer  containing 
many  individual  cells  arranged  in  an  appropriately  symmetric  geometry. 

As  a charged  particle  passes  through  the  chamber,  it  interacts  electromagnet- 
ically  with  the  gas  molecules  in  the  chamber  in  a process  called  ionization;  that  is, 
an  incident  charged  particles  ionizes  gas  molecules  along  its  trajectory.  The  electrons 
from  the  ionized  gas  molecules  then  “drift”  in  the  electric  field  toward  the  nearest 
anode  wire.  As  the  electron  gets  very  close  to  the  anode,  the  electric  field  becomes 
very  strong  which  causes  an  avalanche  whereby  further  ionization  is  induced.  The 
result  of  the  avalanche  is  a tremendous  number  of  electrons  collapsing  upon  the  wire 


^ Anode  and  cathode  wires  are  often  called  “sense”  and  “field”  wires,  respectively.  I will  use  anode 
and  sense,  and  cathode  and  field,  interchangeably. 
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Figure  2.8:  Cross  Section  of  the  Precision  Tracking  Layers. 

in  a very  short  amount  of  time  (less  than  one  nanosecond).  When  this  happens  to  a 
sense  wire,  we  say  that  the  wire  is  “hit”. 

The  current  on  the  anode  wire  from  the  avalanche  is  collected  and  amplified  by 
electronics  at  one  or  both  ends  of  the  anode  wire.  Typically,  an  ADC  and  a TDC 
is  found  among  the  components  of  a drift  chambers’  readout  electronics.  The  ADC 
measurement  corresponds  to  the  amount  of  charge  collected  on  the  wire,  and  the  TDC 
measurement  corresponds  to  the  “time”  of  the  hit.  A calibration  is  determined  that 
maps  the  ADC  values  into  specific  ionization  measurements  of  the  incident  particle, 
thereby  providing  some  information  on  the  particle’s  mass  and  species.  A calibration 
is  also  found  to  turn  the  TDC  values  into  drift  distances;  that  is,  the  distance  of 
closest  approach  an  incident  particle  makes  to  the  anode  wire  (see  [16]  for  further 
discussion). 

The  innermost  tracking  chamber  in  CLEO  II  is  the  Precision  Tracking  Layers, 
also  known  as  the  “PT”  or  “PTL”.  A cross  section  of  the  PTL  is  shown  in  Fig.  2.8. 
The  PTL  is  a six  layer  device,  with  64  straw  tubes  per  layer.  Each  tube  has  an  anode 
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Figure  2.9:  The  Vertex  Detector  anodes. 
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wire  going  down  the  middle  of  it.  The  anode  wires  are  gold-plated  tungsten,  with 
a diameter  of  15  fj,m.  The  tubes  themselves  are  compose  of  aluminized  mylar.  The 
tubes  are  all  glued  together,  with  pieces  of  GlO  plastic  functioning  as  endplates  to 
hold  the  anodes.  During  CLEO  II,  there  were  two  gases  used  in  the  PTL  tubes.  At 
first,  a 50-50  mix  of  argon  and  ethane  was  used.  In  1992,  the  gas  was  switched  to 
dimethyl  ether,  which  gives  better  position  resolution.  The  PTL  sense  wires  are  only 
capable  of  making  r — (j)  measurements. 

After  exiting  the  PTL,  the  next  device  a particle  encounters  is  the  Vertex  De- 
tector (VD).  The  layout  of  the  VD  wires  are  shown  (with  the  PTL  inside)  in  Fig.  2.9. 
The  VD  is  a bit  different  than  the  PTL  - it  is  not  a straw  tube  device.  Wires  are  used 
instead  of  tubes  to  define  the  boundaries  of  the  drift  cells.  This  reduces  the  amount 
of  material  an  incident  particle  has  to  go  through  (and  is  a lot  cheaper) , but  has  more 
complicated  isochrones  (an  isochrone  is  the  collection  of  points  in  a cell  whose  drift 
times  are  the  same).  There  are  10  layers  of  hexagonal  drift  cells  in  the  VD,  with  64 
cells  per  layer  in  layers  1-5,  and  96  cells  per  layer  in  layers  6-10,  for  a total  of  800  drift 
cells  in  all.  The  sense  wires  are  composed  of  a nickel-chromium  alloy,  and  field  wires 
of  gold-plated  aluminum.  The  diameter  of  the  sense  wires  is  20  ^m.  There  are  two 
thicknesses  of  field  wires  used  in  the  VD;  the  thicker  162.5  /xm  diamter  wire  is  used 
in  layers  1,  2,  6,  and  7,  and  142.2  //m  diameter  wire  is  used  everywhere  else.  The 
gas  is  a 50-50  mix  of  argon-ethane.  Gas  is  introduced  into  the  VD  through  through 
four  “gas  holes”  in  layer  5.  The  cells  occupied  by  the  gas  holes  do  not  contain  sense 
wires.  The  endplates  through  which  the  wires  are  threaded  and  gas  introduced  are 
again  made  of  GlO,  and  the  inner  and  outer  walls  of  the  device  are  carbon  filament 
tubes. 

The  VD  sense  wires  are  instrumented  on  both  ends  of  the  wire.  This  allows 
2 measurements  to  be  made  via  the  “charge  division”  method.  For  further  details 


41 


Outer 

Cathode 

Strips 


Inner 

Cathode 

Strips 


Carbon 

Filament 

Tube 


Figure  2.10:  The  Vertex  Detector  cathodes. 

Table  2.1:  Parameters  of  Vertex  Detector  cathode  layers. 

Layer  # of  (/>  segments  # of  2 segments  Pad  size  in  ^ (mm) 

VD  Inner  8 64  5.85 

VD  Outer  8 96  6.85 

regarding  charge  division  see  reference  [17].  The  VD  sense  wires  of  course  can  make 
r — 4>  measurements,  as  well. 

The  VD  also  has  segmented  cathode  strips  on  the  surface  of  the  inner  and  outer 
walls  of  the  chamber.  These  are  shown  in  Fig.  2.10.  The  segmentation  is  in  (j)  and  z, 
though  it  is  much  finer  in  2:  in  order  to  make  ^ position  measurements  along  a charged 
particle’s  trajectory.  The  cathode  strips  also  serve  to  shape  the  sense  wire’s  field  cage 
in  layers  1 and  10  on  the  inner  and  outer  edges,  respectively.  A table  parameters 
describing  the  VD  cathode  layers  is  displayed  in  Table  2.1. 

The  outermost  tracking  chamber  is  the  main  drift  chamber,  or  just  DR.  The 
layout  of  the  DR  anode  wires  are  shown  in  Fig.  2.11.  It  is  much  large  in  volume 
than  the  inner  two  tracking  devices.  A total  of  12,240  sense  and  36,270  field  wires 
are  arrange  in  51  layers  of  roughly  square  cells.  The  number  of  cells  per  layer  is 


X □ 
X X 
X □ 


Outer  Shell 


□xnxnxDXDxaxDxnxnxDxnxn 
xxxxxxxxxxxxxxxxxxxxxxx 
xnxnxDxDXDxnxDxnxnxaxnx 
xxxxxxxxxxxxxxxxxxxxxxx 
□xDxnxnxnxnxnxnxnxnxnxn 
j^xxxxxxxxxxxxxxxxxxxxxx 
xnxDXDxDxnxDxnxnxDXDxax 
xxxxxxxxxxxxxxxxxxxxxxx 
□xnxaxnxDxnxnxDxnxDxnxn 
vxxxxxxxxxxxxxxxxxxxxxv 
bxbxbxbxbxbxbxbxbxbx ■ X ■ 
xxxxxxxxxxxxxxxxxxxxx 

□xnxDXDxnxDxnxnxnxnxn 

j^xxxxxxxxxxxxxxxxxxxx 

jjQxDxnxnxDxnxnxnxDxnx 

xxxxxxxxxxxxxxxxxxxxx 

□xaxDxnxDXDxnxnxnxnxn 

xxxxxxxxxxxxxxxxxxxxx 

jjQxnxnxnxnxoxnxnxnxnx 

xxxxxxxxxxxxxxxxxxxxx 

□xDxDxnxnxnxnxnxnxnxn 

vxxxxxxxxxxxxxxxxxxxx 

gxaxBXBxaxBXBxaxixaxB 


^ Outer  Cathode 


□ Axial  Sense  Wires 
B Stereo  Sense  Wires 
X Field  Wires 


a X a X 

X X X X 
X □ X □ 
X X X X 

□ X □ X 
X X X X 
X □ X □ 

: X X X X 

X a X a 

< X X X 

□ X □ X 
X X X X 
X □ X □ 
X X X X 
□ X □ X 

: X X X X 
X B X a 
c X X X 
X □ X □ 

X X X X 

□ X □ X 
X X 


XOXDXOX  *Xx 

XXXXXXXxv  * 

X ■ X ■ X . X V X X 

XXXXXxxx  ’ 

^ ^ X ^ 


□ X □ 


□ X □ X _ 

xxxxx^i'^^a 

xnxDxDxnJ' 

□xnxDxn;^x 
XXXXXXxx^^^ 
X B X a X a x V 

xaxnxn^><x 

xxxxxxx^ °x 
xnxDxQ^ *Xx 


Inner  Cathode 
Carbon  Filament 


Vertex  Detector 


Figure  2.11:  The  Drift  Chamber  anodes. 
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Figure  2.12:  The  Drift  Chamber  cathodes. 


such  that  the  cell  size  is  approximately  constant.  Forty  of  the  51  layers  have  their 
sense  wires  parallel  to  the  beam  line  {z  axis),  but  11  layers  are  so-called  “stereo” 
layers  and  have  their  sense  wires  strung  at  a small  “stereo”  angle  to  the  2:  axis.  The 
stereo  layers  are  then  able  to  provide  some  measurement  of  the  z position  along  a 
track’s  path.  The  down  side  of  the  stero  angle  approach  is  that  position  resolution  is 
easily  harmed  by  the  fact  that  the  cell’s  isochrones  are  a function  of  2.  This  makes 
calibration  of  the  drft  time-to-distance  relationship  more  difficult  (and  it  is  already 
difficult).  The  sense  wires  are  20  /xm  diameter  gold-plated  tungsten.  There  are  again 
two  types  of  field  wires.  Field  wires  in  DR  layers  1 through  40  are  110  /xm  diameter 
gold-plated  aluminum,  and  in  layers  41  through  51  are  110  /xm  diameter  gold-plated 
copper-beryllium.  The  wires  are  supported  by  3.175  cm  thick  aluminum  endplates, 
and  the  gas  in  the  chamber  is  again  the  50-50  mix  of  argon-ethane.  The  inner  wall 
of  the  chamber  is  a carbon  filament  reinformed  epoxy  tube,  and  the  outer  wall  is  a 
composite  of  panels. 

As  in  the  VD,  segmented  cathode  surfaces  shape  the  field  cage  on  the  inner 
surface  of  the  innermost  layer  and  on  the  outer  surface  of  the  outermost  layer.  The 


Table  2.2:  Parameters  of  DR  cathode  layers. 

^ of  0 segments  # of  2:  segments  Pad  size  in  z (mm) 
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Layer 


DR  Inner 

16 

96 

10.00 

DR  Outer 

8 

192 

9.50 

DR  cathodes  are  shown  in  Fig.  2.12,  and  a table  describing  their  parameters  is  shown 
in  Table  2.2. 

2. 3. 1.3  Time  of  Flight  System 

The  time  of  flight  system  is  composed  of  a series  of  plastic  scintillators  and 
is  arranged  as  two  endcaps  and  a barrel.  As  the  name  suggests,  the  purpose  of  this 
systme  is  to  measure  the  time  a particle  spends  “in  flight” ; i.e.  the  time  a particle  takes 
to  travel  from  the  interaction  point  to  the  scintillators.  Time  of  flight  information  is 
used  in  several  ways;  to  trigger  the  detector,  to  identify  the  which  car  in  the  train  a 
collision  came  from,  and  to  provide  a means  of  particle  identiflcation. 

In  the  barrel  section  there  are  64  scintillators,  each  10  cm  wide  and  279.4  cm 
long.  The  angular  coverage  of  the  barrel  region  is  9 between  35  and  145  degrees. 
In  each  endcap  there  are  28  pie-slice  shaped  scintillators  arranged  around  the  beam 
pipe  covering  polar  angles  between  18  and  37  degrees  on  the  west  side,  and  143  to 
162  degrees  on  the  east.  The  scintillators  are  all  about  5 centimeters  thick.  The  tiny 
amount  of  overlap  between  the  barrel  and  endcap  systems  can  be  seen  most  easily  in 
Fig.  2.7. 

2. 3. 1.4  Crystal  Calorimeter 

The  crystal  calorimeter  is  composed  of  7800  thallium-doped  Csl  crystals,  and 
like  the  time  of  flight,  is  arranged  as  two  endcaps  and  a barrel.  The  barrel  contains 
6144  crystals,  and  each  one  is  pointed  towards  the  interaction  point  so  that  photons 
originating  there  strike  the  barrel  crystal  faces  at  nearly  normal  incidence  (this  can 
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Figure  2.13:  Crystal  Calorimeter  Barrel. 


be  seen  in  Figures  2.6  and  2.7).  The  828  crystals  in  each  endcap  are  aligned  with  the 
z axis.  Each  crystal  is  approximately  5 cm  square  and  30  cm  long. 

Photodiodes  are  attached  to  the  back  of  each  crystal  to  convert  light  into  elec- 
trical signals.  For  redundancy,  there  are  four  photodiodes  per  crystal.  When  the 
crystals  are  read  out,  a mixer/shaper  card  sums  the  four  signals  from  each  crystal 
and  shapes  the  signal  for  input  into  an  ADC. 

A section  of  the  barrel  crystal  holder  is  shown  in  Fig.  2.13.  A quadrant  of  an 
endcap  crystal  holder  shown  in  Fig.  2.14. 


2.3. 1.5  Muon  System 

Like  the  crystal  calorimeter  and  time  of  flight  systems,  the  muon  system  is 
arranged  as  two  endcaps  and  a barrel.  The  purpose  of  the  muon  system,  as  the  name 
implies,  is  to  identify  the  passage  of  a muon. 

The  muon  detectors  are  plastic  streamer  counters  that  are  approximately  5 m 
long  and  8.3  cm  wide.  The  counters  are  placed  in  between  layers  of  iron,  which  only 
mouns  can  penetrate  with  high  efficiency.  In  each  gap  between  the  iron,  there  are 
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Figure  2.14:  Crystal  Calorimeter  Endcap. 


three  layers  of  counters.  Each  counter  contains  eight  compartments  containing  an 
anode  wide  and  is  filled  with  a 50-50  mix  of  argon-ethane  gas.  The  anode  wires  are 
50  /rm  in  diameter,  and  are  made  of  silver-plated  copper-beryllium.  Each  anode  is 
surronded  on  three  sides  with  graphite  cathodes,  and  within  a counter,  all  the  anode 
wires  are  ganged  together.  When  a hit  is  recorded,  the  anode  wire  position  provides 
the  ({)  coordinate  of  the  hit,  and  charge  division  is  used  to  extract  the  2;  coordinate. 

2.3.2  The  CLEO  II.V  Detector 

After  a several  month  shutdown  for  the  installation  of  the  Phase  II  CESR/CLEO 
Upgrade,  CLEO  II.V  began  taking  data  in  November  1995.  This  upgrade  included 
the  replacement  of  the  beam  pipe  as  mentioned  in  Section  2. 1.3.1,  as  well  as  the 
removal  of  the  PT  and  installation  of  a new  tracking  device  to  take  its  place:  the 
Silicon  Vertex  Detector  (SVX).  The  SVX  is  described  in  reference  [19].  A diagram  of 
it  is  shown  in  Fig.  2.15. 

The  SVX  is  not  a drift  chamber;  it  is  a silicon  strip  detector,  and  the  means 
of  operation  are  quite  different  from  those  of  the  tracking  devices  we  have  discussed 
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Figure  2.15:  The  CLEO  II. V Silicon  Vertex  Detector. 
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so  far.  The  SVX  is  a three  layer  device  consisting  of  96  silicon  wafers,  each  300  iim 
thick.  As  shown  in  Fig.  2.15,  the  wafers  are  arranged  into  octants;  each  octant  has 
12  wafers.  The  octants  can  be  seen  in  the  r — cj)  view  to  overlap  each  other  slightly 
in  the  corners,  and  results  in  tracks  which  can  have  hits  on  two  different  wafers  in  a 
layer.  This  is  useful  in  aligning  the  detector.  A thorough  summary  of  the  parameters 
describing  the  SVX  can  be  seen  in  Table  2.3. 

The  wafers  are  instrumented  and  read  out  on  both  sides;  hence  the  name  “Dou- 
ble Sided  Silicon  Detector”  (DSSD)  used  in  Fig.  2.15.  The  instrumentation  on  each 
side  consists  of  an  array  of  aluminum  strips  on  the  wafer  surface.  These  strips  are 
ultimately  connected  to  preamplifiers  stationed  at  the  end  of  the  detector.  On  the 
top  of  the  wafer  (facing  away  from  the  beam  line),  the  aluminum  strips  sit  atop  p 
type  implants,  with  strips  and  implants  parallel  to  the  beam.  On  the  bottom  of  a 
wafer  (facing  the  beam  line),  the  strips  connect  to  n type  implants,  with  strips  and 
implants  perpendicular  to  the  beam.  The  entire  wafer  forms  a,  p — n junction.  When 
reverse  bias  is  applied  across  the  wafer,  a sensitive  region  depleted  of  mobile  charge 
is  formed. 

Charged  particles  traversing  the  wafer  lose  energy  as  in  any  other  material. 
In  the  sensitive  region  of  a silicon  wafer,  this  lost  energy  is  used  to  create  electron- 
hole  pairs.  Approximately  3.6  eV  is  required  to  create  a single  electron-hole  pair, 
much  lower  than  the  ionization  energy  in  drift  chamber  gases  (40  to  100  eV).  For 
minimum  ionizing  particles  in  the  CLEO  SVX,  this  translates  to  about  21,000  pairs. 
The  electrons  and  holes  then  travel  (in  opposite  directions)  in  the  electric  field  applied 
by  the  bias  to  the  surfaces  of  the  wafers  and  end  up  on  the  strips.  In  this  way,  the 
appearance  of  charge  on  a strip  (a  “hit”  strip)  indicates  the  passage  of  a charged 
particle.  Since  the  strips  on  the  top  of  the  wafer  {p  side)  are  parallel  to  the  beam, 
they  provide  measurements  of  the  (j)  location  of  the  track  at  the  center  of  the  wafer 
(this  is  why  the  top  of  the  wafer  is  called  the  r — 4>  side) . Similarly,  the  strips  on  the 
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bottom  of  the  wafer  [n  side)  measure  the  z position  and  motivate  the  term  z side 
referring  to  the  bottom  of  a wafer. 

Previously,  we  mentioned  that  the  readout  electronics  are  mounted  on  the  ends 
of  the  detector.  On  the  r — 0 side,  the  strips  run  parallel  to  the  beam  line,  mak- 
ing readout  a straightforward  process.  However,  the  z side  strips  themselves  are 
perpendicular  to  the  beam.  Reading  them  out  directly  would  require  the  front  end 
electronics  to  be  placed  in  the  tracking  volume,  leading  to  increased  multiple  scat- 
tering as  charged  particles  interacted  with  the  material  of  the  electronics.  Increased 
multiple  scattering  results  in  degradation  of  tracking  resolution.  So  instead,  the  2: 
side  strips  are  connected  to  traces  which  carry  the  signals  to  the  ends  of  the  detector. 
Each  trace  is  connected  three  strips  thereby  “ganging”  the  strips  together.  Gang- 
ing leads  to  ambiguity  in  which  strip  actually  collected  charge,  and  therefore  makes 
tracking  more  difficult.  The  ganged  strips  are  spaced  far  apart  to  minimize  this  effect. 
However,  this  significantly  reduces  the  number  of  channels  the  front  end  electronics 
has  to  deal  with. 

In  all,  the  SVX  contains  26,208  “read  out”  strips,  or  active  channels.  This  is 
more  than  the  rest  of  the  channels  in  the  entire  detector  combined.  To  read  out 
the  entire  detector  every  event  would  slow  down  the  data  aqcuisition  system  to  a 
crawl.  Therefore,  the  events  are  sparsified  online  which  drops  channels  which  do  not 
have  any  charge.  Unlike  the  sparsifiers  used  in  the  other  detector  subsystems  which 
are  implemented  in  hardware,  the  SVX  sparsifier  is  a program  that  runs  on  the  same 
Motorola  68040  CPU  chips  found  in  older  Apple  and  Amiga  personal  computers.  The 
advantage  of  running  a software  sparsifier  is  that  its  algorithm  can  be  changed  easily. 
The  disadvantage  is  that  it  is  not  as  fast  as  a sparsifier  implemented  in  hardware. 

Other  changes  resulting  the  Phase  II  upgrade  included  changing  the  gas  in 
the  DR  from  the  50-50  mix  of  argon-ethane  to  60-40  helium-propane.  The  helium- 
propane  mixture  has  better  position  resolution  due  to  the  greater  number  of  primary 
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ions  created  per  unit  length  of  cell  traversed  (bigger  hydrocarbon  molecules).  In 
addition,  hit  efficiency  at  the  edge  of  a drift  cell  was  improved  [20].  Another  change 
was  the  installation  of  multi-hit  electronics  on  one  side  of  the  VD  as  part  of  research 
and  development  for  the  Phase  III  upgrade  [21].  Though  these  electronics  were  not 
used  for  tracking,  their  presence  (or  rather  the  removal  of  the  old  electronics)  did 
have  some  effect  on  CLEO  II.  V tracking  performance;  the  multihit  electronics  do  not 
have  ADCs.  This  makes  the  method  of  charge  division  impossible  for  assigning  a 2 
measurement  to  a VD  hit.  In  CLEO  II,  tracking  of  low  momentum  charged  particles 
often  relied  on  charge  division  hits  to  provide  this  z information. 

The  major  effect  in  the  Phase  II  upgrade,  however,  was  brought  on  by  installa- 
tion of  the  SVX.  Exchanging  the  PT  for  the  SVX  provided  significant  improvements 
to  the  tracking  system.  This  is  for  two  reasons.  Eirst,  the  SVX  measures  both  the 
r — cf)  and  the  2 position  of  a track  with  a single  layer.  The  PT,  on  the  other  hand,  had 
no  ability  to  make  a 2:  measurement  at  all.  Second,  the  position  resolution  afforded 
by  the  SVX  is  much  better  than  the  PT,  even  mr  — (j)  where  the  PT  was  already  the 
highest  resolution  tracking  device,  far  better  than  the  VD  or  DR.  The  improvement 
in  the  z and  0 resolution  on  a charged  particle’s  trajectory  is  dramatic.  The  result  is 
that  after  the  Phase  II  upgrade,  CLEO  is  in  an  excellent  position  to  make  precision 
lifetime  measurements  short-lived  particles. 
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Table  2.3:  Silicon  Vertex  Detector  parameters. 


Layer  1 

Layer  2 

Layer  3 

Radius  to  inner  surface 

23.5  mm 

32.5  mm 

46.9,48.1  mm 

Total  length 

65.580  mm 

81.270  mm 

2x65.580  mm 

Total  width 

22.568  mm 

30.128  mm 

2x22.568  mm 

Active  length 

59.535  mm 

76.356  mm 

2x59.525  mm 

Active  width 

21.168  mm 

28.890  mm 

2x21.168  mm 

cos^  coverage 

0.92 

0.92 

0.92 

Thickness 

300  /um 

300  fim. 

300  /xm 

Number  of  DSSD  detectors 

16 

16 

64 

Number  of  CAMEX  chips 

96 

128 

96 

Number  of  JAMEX  chips 

0 

0 

96 

Readout  bond  pad  pitch 

97  //m 

97  /xm 

97  /xm 

Number  of  active  channels 

6048 

8064 

12096 

Detector  Outer  Side  (r-0, 

a = 12/xm) 

Implant 

p-type 

Number  of  implants 

753 

1005 

3012 

Implant  width 

10  iJ,m 

10  /xm 

10  /xm 

Implant  pitch 

28  fim 

28.5  /xm 

28  /xm 

Readout  width 

5 /um 

5 /xm 

5 /xm 

Readout  pitch 

112  //m 

114  /xm 

112  /xm 

Readout  strips 

189 

252 

378 

Detector  Inner  Side  {r-z, 

a = 30/xm) 

Implant 

n-type 

Number  of  implants 

567 

756 

1134 

Implant  width 

20  jim 

20  /xm 

20  /xm 

Implant  pitch 

105  /xm 

101  /xm 

105  /xm 

Approximate  Barrier  Width 

69  /xm 

65  /xm 

69  /xm 

Readout  width 

5 /xm 

5 /xm 

5 /xm 

Readout  pitch 

105  ^m 

101  /xm 

105  /xm 

Readout  strips 

189 

252 

378 

CHAPTER  3 
ANALYSIS  TECHNIQUE 


In  order  to  measure  the  mixing  parameter  y in  the  system,  we  must 

obtain  a large  statistics  sample  of  and  decays  into  CP  eigenstate  final  states, 
and  measure  a proper  time  for  each  or  parent.  In  this  chapter,  we  describe 
the  data  samples  used,  the  global  event  requirements  imposed,  and  the  selection 
criteria  and  technique  used  to  choose  and  reconstruct  such  a sample.  As  the  purity 
of  the  sample  is  an  important  issue,  it  is  worth  noting  that  a reconstructed  or 
candidate  is  just  a candidate.  We  have  no  way  of  knowing  whether  or  not  it  is 
real.  Therefore,  the  efficiency  of  our  selection  criteria  and  their  ability  to  reject  false 
candidates  (“background”)  is  discussed. 

3.1  Data  Samples  and  Event  Selection 

The  data  used  in  this  analysis  were  collected  by  the  CLEO  II. V detector.  We 
use  full  sample  of  9.0  fb~^  taken  on  and  near  the  T(4S). 

The  data  in  this  sample  is  logically  divided  into  datasets,  each  of  which  corre- 
sponds to  several  months  of  data-taking;  typically  a few  thousand  runs.  A boundary 
between  datasets  is  created  when  changes  that  can  affect  either  detector  acceptance 
or  resolution  are  imposed  (such  as  opening  the  detector  pole  tips  for  maintenance). 

The  most  detailed  monte  carlo  simulation  used  in  this  analysis  is  based  on 
GEANT  [22],  which  incorporates  a full  detector  simulation  and  response.  The  statis- 
tical advantage  over  the  data  is  « 8.25  for  the  KK  and  tttt  samples,  and  ~ 6 for  Kir. 
The  simulation  included  a proper  luminosity-weighted  sampling  of  the  entire  CLEO 
II. V dataset.  This  will  account  for  run-dependent  effects,  such  as  variation  in  beam 
energy  and  calibration  constants,  to  the  best  of  our  ability. 
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Most  of  the  e'^e“  collisions  at  this  energy  do  not  result  in  events  containing 
hadrons,  and  are  not  interesting  to  this  analysis.  Therefore,  we  impose  the  following 
requirements  to  choose  events  likely  to  contain  hadrons: 

• The  number  of  reconstructed  charged  tracks  required  to  be  > 3. 

• If  the  number  of  charged  tracks  is  3 or  4,  the  total  energy  deposited  in  the  calori- 
meter must  be  less  than  65%  of  the  beams’  center  of  mass  energy,  Ecm,  or  the 
highest  energy  shower  in  the  calorimeter  must  be  less  that  60%  of  Et,eam- 

• The  total  visible  energy  must  be  greater  than  20%  of  Ecm- 

• The  total  calorimeter  energy  must  be  greater  than  15%  of  Ecm  or  number  of 
charged  tracks  greater  than  4. 

• Primary  vertex  position  is  required  to  be  within  2 cm  of  origin  in  CLEO  world 
coordinate  system  in  xy,  and  within  7.5  cm  in  z. 

The  above  criteria  retain  over  95%  of  udsc  and  more  than  99%  of  BB  events,  while 
removing  more  than  99%  of  the  ubiquitous  e'*'e“  — > e“*"e“  and  e''‘e~  e^e~'y  events. 

As  the  performance  of  the  tracking  system  is  paramount  in  this  analysis,  we 
demand  that  the  confidence  level  for  the  indentification  of  which  car  of  the  CESR 
bunch  train  the  annihilation  orinated  from  is  more  than  0.5.  This  cut  is  ~ 
100%  efficient  for  hadronic  events.  If  the  wrong  car  was  selected,  this  would  insert  an 
incorrect  offset  of  order  of  tens  of  nanoseconds  (how  much  depends  upon  what  car  the 
event  originated  from  and  what  car  was  mistakenly  identified)  into  the  determination 
of  drift  times,  resulting  in  extremely  poor  tracking  resolution. 

3.2  Building  Blocks 

Excepting  for  some  “uninteresting”  types  of  events,  we  do  not  normally  observe 
(detect)  all  the  particles  created  by  an  annihilation.  Indeed,  these  “uninterest- 
ing” types  of  events  are  what  our  event  selection  criteria  (described  in  Section  3.1) 
are  designed  to  throw  away!  In  the  events  which  survive  our  selection  criteria,  the 
result  of  the  e‘*'e“  annihilation  is  usually  the  creation  of  many  particles.  Some  of  these 
particles  leave  evidence  of  their  existence  inside  CLEO’s  various  detector  elements. 


and  some  do  not. 
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A stable  particle  may  not  be  directly  observable  by  CLEO  for  many  reasons. 
Maybe  the  particle  simply  does  not  interact  sufficiently  with  the  CLEO  detector 
material  in  order  for  its  trajectory  to  be  detected  - this  is  case  with  neutrinos,  for 
example.  Or  perhaps  the  direction  of  the  particle  is  such  that  it  never  traverses  any 
of  CLEO’s  detector  elements  - i.e.,  the  particle  escapes  down  the  beam  pipe.  These 
are  “missing”  particles,  and  it  is  crucial  for  certain  analyses  (not  this  one)  to  employ 
techniques  to  try  to  recover  such  cases  [28].  But  there  is  another  important  reason  a 
particle  may  not  be  directly  detected  by  CLEO  - it  is  not  stable,  and  decays  before 
encountering  any  of  CLEO’s  detectors.  For  this  class  of  particles,  it  may  still  be 
possible  to  “reconstruct”  the  momentum  of  the  unstable  particle  and  position  of  the 
decay  point  - if  the  daughter  particles  are  detected. 

In  a High  Energy  Physics  analysis,  one  build  high  level  objects  from  more  prim- 
itive ones.  An  example  of  this  is  the  construction  of  tracks  from  the  the  low  level 
hit  information  provided  by  the  tracking  chambers.  The  reconstruction  of  particle 
candidates  is  logically  very  similar  - just  as  a track  can  be  thought  of  as  the  “sum” 
of  the  hits  used  to  form  it,  a particle  which  decays  into  a number  of  daughters  can 
be  viewed  as  the  sum  of  those  daughters.  Applying  the  reconstruction  methodology 
iteratively,  complex  decay  chains  and  even  entire  events  can  be  fully  reconstructed. 

In  our  reconstruction,  the  detected  daughters  of  the  decay  are  our  “build- 
ing blocks”  for  the  D^.  For  the  decays  we  are  interested  in,  those  building  blocks 
are  K and  tt  charged  particles.  Their  selection  criteria  is  described  below. 

3.3  Charged  Particle  Selection 

The  CLEO  II. V detector’s  tracking  chambers  generally  provide  excellent  mo- 
mentum and  position  resolution  for  charged  particle  candidates.  However,  measuring 
the  lifetime  of  a particle  so  short-lived  as  the  is  a hard  business,  and  is  near 
the  limit  of  CLEO’s  resolution.  We  therefore  use  high  quality  track  candidates  to 
reconstruct  the  . 
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To  choose  charged  particle  candidates  of  good  quality,  we  require  candidates 
to  be  successfully  Kalman  fit  - this  requirement  ensures  that  the  track  parameters 
have  proper  ionization  energy-loss  corrections  and  the  errors  on  the  track’s  trajectory 
are  reasonable.  Tracks  also  must  have  mean  residuals  in  the  drift  chamber  of  less 
than  1 mm.  As  the  position  resolution  of  a drift  chamber  hit  is  about  180  //m,  this 
requirement  removes  only  very  poorly  measured  tracks. 

As  position  information  is  paramount  in  lifetime  determination,  we  demand 
each  daughter  candidate  have  two  hits  in  each  view  (r0  and  rz)  of  the  highly 
position-sensitive  SVX. 

We  also  do  not  consider  tracks  whose  momentum  is  less  than  250  MeV  as 
daughter  candidates,  as  our  simulation  tells  us  that  there  are  very  few  such  candi- 
dates. Since  low  momentum  charged  particles  tend  not  to  have  their  trajectories 
reconstructed  as  well  as  high  momentum  particles,  this  criteria  also  removes  some 
poorly  measured  candidates. 

“Dredge”  seeds  are  track  candidates  found  by  the  pattern  recognition  in  a last- 
ditch  attempt  to  use  all  hits  in  the  tracking  system.  Any  hits  left  over  from  the  more 
reliable  pattern  recognition  methods  are  “dredged”  into  track  candidates.  These 
candidates  tend  to  be  complete  figments  of  the  tracker’s  imagination.  Tracks  whose 
seeds  are  from  the  “dredge”  level  of  the  pattern  recognition  are  not  considered,  as 
“dredge” -ing  is  prone  to  making  tracks  which  do  not  correspond  to  a real  particle. 
Further,  we  demand  that  the  tracks  be  approved  by  the  track  quality  package  TMNG, 
which  removes  yet  more  pattern  recognition  anomolies. 

Finally,  we  demand  that  charged  particle  candidates  are  not  part  of  a well 
identified  A,  or  7 conversion. 

3.4  Neutral  D Meson  Reconstruction 

A candidate  is  formed  by  combining  the  4-momenta  of  a set  of  vertex  con- 
strained charged  particle  candidates,  which  have  passed  the  selection  criteria  of  Sec- 
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tion  3.3;  we  consider  all  possible  pairs  of  oppositely  charged  tracks  under  the  7r'''7r“, 
K^'k~ , , and  K~K^  hypotheses.  A covariance  matrix  describing  the  errors 

on  the  momentum  and  decay  vertex  position  is  calculated  from  the  covariance 
matrices  of  the  daughters. 

To  reduce  the  number  of  candidates  from  random  combinations  of  tracks, 
we  seek  out  candidates  which  come  from  the  decay  D*^  — >•  . To  build  the 

D*  candidate,  we  trace  the  direction  of  the  candidate’s  trajectory  back  from  the 
reconstructed  decay  vertex  to  the  beam  spot.  The  {x,y,z)  coordinates  along  the 
trajectory  that  are  most  consistent  with  the  beam  (using  the  beam’s  covariance 
matrix)  are  identified,  and  a tTs  candidate  is  forced  through  this  point.  The  selection  of 
the  TTs  candidates  is  the  same  as  outlined  in  Section  3.3,  except  that  the  momentum 
requirement  is  lowered  to  50  MeV,  and  we  do  not  demand  any  SVX  information. 
A D*  candidate  is  the  built  from  the  refit  tTj  and  the  undisturbed  trajectories. 
This  procedure  improves  the  resolution  on  the  mass  difference.  Am  = m£>-  — tud 
dramatically.  We  demand  that  |Am|  = \mo*  — tud]  is  within  1 MeV  of  it’s  nominal 
value  of  145.4  MeV. 


A cartoon  depicting  of  the  construction  of  and  D*  candidates  is  shown  in 
Figure  3.1. 

In  order  to  ensure  candidates  of  good  quality  and  to  reduce  background 
from  random  charged  particles,  we  require  the  probability  of  vertex  chi-squared  for 
the  decay  vertex  be  greater  than  10“^. 

In  CLEO’s  environment,  a flies  on  order  of  200  /xm  before  decaying.  This 
is  not  very  far,  and  means  that  the  decay  vertex  is  very  near  to  the  high  number 
of  tracks  emerging  from  the  event’s  primary  vertex  (where  the  e'*'e“  collision  took 
place).  This  means  that  many  charged  particles  originating  from  the  primary  vertex 
pass  our  vertex  quality  selection  criteria  whether  they  come  from  from  a D meson  or 
not.  This  kind  of  combinatoric  background  is  mostly  populated  by  low  momentum 
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Silicon  Detector 


Reconstructed  track 
without  refit 


Beam  Spot 


^ Slow  Pion 


D decay  vertex 


Primary  Vertex 


Other  charmed 
hadron 

Figure  3.1:  Diagram  of  D and  D*  candidate  reconstruction, 


particles  (whose  errors  are  largest,  and  are  therefore  most  likely  to  be  “happy”  to  be 
constrained  to  the  slightly  displaced  D meson  vertex).  We  can  remove  most  of  this 
background  by  cutting  on  |cos(0*)|,  where  cos{9*)  is  defined  by: 


cos (9*)  = 


Pdau  ' PD 


\Pdau\\P 


D 


(3.1) 


and  is  the  momentum  of  one  of  the  daughter  tracks  (we  choose  the  one  whose 
charge  matches  that  of  the  candidate)  in  the  rest  frame  of  the  D candidate.  That 
is,  9*  is  the  angle  between  the  D direction  in  the  lab  and  the  direction  of  the  chosen 
daughter  in  the  rest  frame  of  the  D.  We  expect  that  the  distribution  of  |cos(0*)|  for  the 
signal  is  fiat,  and  that  for  the  background  is  peaked  near  one  with  the  peak  populated 
by  random  low  momentum  tracks  which  happen  to  be  moving  along  (or  against)  the 
direction  of  the  D candidate’s  boost.  This  is  indeed  what  we  find.  Demanding  that 


cos(0*)|  < 0.8  significantly  reduces  this  background. 
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Recall  that  CLEO/ CESR  spends  approximately  two-thirds  of  its  time  running 
on  the  T(4S')  resonance,  which  can  decay  to  a pair  of  BB  mesons.  Though  the 
T(45)  is  just  above  threshold  for  BB  production,  the  resulting  B mesons  are  not 
produced  at  rest,  and  on  average  fly  about  30  microns  before  decaying.  While  30 
microns  is  quite  a short  distance,  it  is  not  completely  negligible  on  the  scale  of  the 
flight  distance  error  on  a reconstructed  As  a determination  of  the  candidate 
lifetime  requires  precise  knowledge  of  both  the  production  and  decay  vertices,  we 
therefore  wish  to  exclude  D“s  originating  from  a B meson  decay.  At  CLEO,  D°s 
originating  from  B meson  decay  (directly  or  indirectly)  have  a strikingly  different 
momentum  distribution  than  D^s  produced  by  other  means.  This  provides  us  with  a 
powerful  means  to  discriminate  against  candidates  from  B mesons  (see  Figure  3.2); 
we  require  that  the  momentum  of  the  candidate  be  greater  than  2.3  GeV/c.  This 
criteria  removes  candidates  with  poorest  proper  time  resolution,  so  the  significant 
loss  of  efficiency  does  not  significantly  affect  the  resolution  on  the  fitted  lifetimes. 

Finally,  we  demand  that  for  each  candidate  hypothesis  (tt+tt", 

and  K~K'^),  we  demand  the  invariant  mass  of  all  the  other  hypotheses  is 
more  than  4 cr  inconsistent  with  the  nominal  D mass.  This  primarily  keeps  sym- 
metric Ktt  decays  out  of  our  sample  (i.e.,  A'+tt"  and  JT“7r+  candidates  formed  with 
the  same  two  tracks),  where  the  assignment  of  tracking  hypothesis  kinematically  is 
not  obvious.  It  also  keeps  all  our  samples  independent;  that  is,  no  pair  of  D decay 
hypotheses  uses  the  same  distinct  tracks. 

3.5  Proper  Time  Measurement 

The  proper  time  elapsed  between  the  production  and  decay  of  a neutral  particle 
is  given  by 

, ^ {T dec  ^prod^  ’ P 

C 1^ 


(3.2) 


Fraction  of  Particles  / 50  MeV/c 
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Figure  3.2:  Neutral  D meson  momentum  spectrum  from  CLEO  II. V monte 

carlo  for  D mesons  originating  from  B meson  events,  and  for  D 
mesons  originating  from  udsc  events  (continuum). 
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where  m is  the  particle’s  mass,  f^ec  is  it’s  decay  vertex  position,  Tprod  is  it’s  production 
point,  and  p is  it’s  momentum  vector.  We  have  m,  fdec,  and  p from  our  D meson 
reconstruction.  What  is  left  to  estimate  is  the  production  vertex,  fprod- 

The  beam  position  is  an  obvious  candidate  for  fprod,  as  is  the  result  of  some 
primary  vertex  finding  algorithm,  perhaps  using  the  beam  position  and  size  as  a 
constraint.  We  elect  not  to  use  a primary  vertex  finding  algorithm  due  to  the  in- 
troduction of  bias  in  the  D decay  length  created  when  a track  not  belonging  to  the 
primary  vertex  is  used  to  calculate  it.  If  we  build  a primary  vertex  which  contains 
any  of  the  daughter  tracks  of  a charmed  meson  on  the  “other  side”  of  the  event,  the 
calculation  of  the  vertex  would  be  pulled  away  from  its  true  value.  The  resulting 
bias  is  systematic,  consistently  farther  away  from  the  D meson  decay  vertex  than  the 
truth,  and  leads  to  corresponding  systematic  biases  in  measured  proper  time.  It  is 
difficult  in  practice  to  avoid  this  kind  of  mistake,  due  to  the  small  distances  travelled 
by  D mesons  and  our  finite  detector  resolution.  The  beam  position,  on  the  other 
hand,  is  an  unbiased  (if  not  optimal)  estimator  of  an  event’s  origin.  Further,  the 
beam  position  gives  excellent  resolution  on  the  y coordinate  of  fprod  (~  10  pm)  due 
to  its  very  small  size  in  that  dimension,  and  the  beam  position  is  uncorrelated  with 
our  D meson  candidate.  So  we  elect  to  use  the  beam  position  as  our  estimator  of 

^prod' 

Rather  than  bend  the  reconstructed  D meson  trajectory  through  fprod,  we  in- 
stead follow  the  D meson  trajectory  to  each  of  the  coordinates  of  fprod  separately.  We 
then  have  three  measurements  of  the  proper  time: 
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and  their  associated  errors: 
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where  we  have  used  the  fact  that  fprod  is  uncorrelated  with  our  D meson  candidate. 

These  three  correlated  proper  time  measurements  can  be  combined  into  a single 
optimal  estimator  of  the  D meson  candidate’s  proper  time  and  uncertainty  [23] 

^ ^ WiU 

i 

(3.5) 

ij 

where  the  indices  i and  j run  over  the  coordinates,  and 

* 

{Vt)ij  = {StiStj)  (3.6) 

are  the  weights  and  covariance  matrix  elements,  and 

Sti  = ti  ( ^ - ^ + ^fd££iZ^fPT£f!i  ) (3.7) 
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In  this  analysis,  we  use  t and  aj  as  our  proper  time  and  proper  time  uncertainty, 
respectively.  Note  that  because  of  our  choice  of  fprod,  the  contribution  to  t from  the  2: 
coordinate  is  negligible  due  to  the  large  size  of  the  beam  in  this  coordinate  (~  3 cm; 
this  would  not  be  true  if  we  had  better  knowledge  of  Zprod)'-,  we  include  it  anyway. 

3.6  Lifetime  Measurement 

We  extract  the  lifetime  of  our  D meson  signals  with  an  unbinned  maximum 
likelihood  fit.  The  inputs  to  the  fit  are  the  measured  proper  time,  U,  and  proper  time 
uncertainty,  and  the  reconstructed  invariant  mass,  rrii.  The  mass  is  converted  to 
a signal  probability,  psig^  > which  is  determined  by  a one  dimensional  binned  fit  to  the 
full  invariant  mass  distribution  of  the  input  sample. 

In  the  proper  time  fit,  the  full  likelihood  function  which  is  maximized  has  the 
following  form: 


— -^Sig  + ■^BG 


That  is,  the  total  likelihood  is  the  sum  of  two  separable  pieces,  one  for  signal  and  one 
for  background.  Similarly,  the  background  is  divided  into  two  components: 

= -^BGo  + ■^BGr 

one  component  with  lifetime  and  one  without.  Generally  speaking,  the  components 
of  the  likelihood  containing  lifetime,  £sig  ^-nd  £BGr>  ^re  modelled  as  exponentials 
convoluted  by  a resolution  function.  The  zero  lifetime  background  piece,  £bGo>  is  a 
delta  function  smeared  by  a resolution  function. 

The  signal  likelihood  resolution  function  is  the  sum  of  three  gaussians.  The 
first  gaussian  is  one  whose  width  is  the  measured  proper  time  error,  at^,  times  some 
scale  factor,  5sig  ; the  5sig  is  meant  to  account  for  a possible  uniform  “mistake”  in  the 
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covariance  matrix  elements  of  the  D meson  (and  its  daughters),  as  would  be  the  case  if 
there  was  an  imperfect  description  of  the  detectors  material  used  during  track  fitting. 
For  the  other  two  gaussians,  the  measured  proper  time  errors  are  simply  ignored  and 
the  widths,  crsigi  and  crsig2,  and  their  normalizations,  /sigj  and  /sig2>  can  be  fit  for 
directly.  In  practice,  cTsig2  is  fixed  to  “something  large”;  we  typically  use  8 ps.  This 
helps  convergence.  The  purpose  of  these  two  gaussians  is  to  model  the  proper  time 
smearing  when  at,i  does  not  represent  the  spread  of  the  proper  time  distribution,  as 
would  be  the  case  for  hard  multiple  scattering  or  tracking  system  pattern  recognition 
error  of  one  or  more  of  the  D meson  daughters  (we  refer  to  these  components  as 
“mis-measured”).  The  sum  of  these  three  components  to  the  likelihood  is  normalized 
by  the  total  signal  probability,  Psigp  as  determined  by  the  fit  to  the  invariant  mass 
distribution.  We  further  note  that  if  we  understand  our  detector  well,  we  expect  to 
find  that  the  scale  factor  used  in  the  first  gaussian  is  close  to  1 and  the  fraction  of 
the  signal  in  the  other  gaussians  is  near  zero. 

The  treatment  of  the  background  is  similar  to  that  of  the  signal.  The  total 
background  likelihood  is  normalized  by  the  background  probability,  which  is  just 
(1  — Psigi)-  The  likelihood  considers  two  types  of  background:  background  with  some 
lifetime  tbg,  and  background  with  zero  lifetime.  We  constrain  both  kinds  of  back- 
ground to  have  the  same  resolution  function.  The  model  for  the  resolution  function 
is  the  same  as  for  the  signal  (three  gaussians).  The  first  gaussian,  whose  width  is 
times  some  scale  factor,  5bG)  models  the  “correctly  reconstructed”  background, 
and  the  other  two  whose  width  and  norms  we  fit  for  describe  the  “misreconstructed” 
background.  In  some  sense,  most  background  is  “misreconstructed”  - it  is,  after  all, 
background.  We  will  assign  a systematic  uncertainty  due  to  our  background  model 
later. 

There  are  a lot  of  fit  parameters  here;  rsig  is  obviously  the  most  important.  It 
should  also  be  noted  that  in  one  should  not  read  too  much  into  the  meaning  of  the 
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parameter  tbG)  as  the  component  of  the  background  it  describes  generally  contains 
more  than  one  lifetime.  The  full  likelihood  function  is: 
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where  G{t\a)  = l/{V^a)  exTp{—t^/2a^)  and  E{t\T)  = 1/r  exp(— t/r).  The  meanings 
for  the  fit  parameters  and  inputs  from  the  data  are  described  in  Table  3.1. 

The  lifetime  extraction  procedure  described  here  is  implemented  in  a software 
package  called  LifeFit  [24]. 

We  note  for  the  record  that  the  parameterization  of  the  likelihood  function  is 
different  from  and  more  general  than  used  in  the  published  CLEO  II.V  measurement 
of  the  Kn  lifetime  [25].  The  difference  is  that  in  the  Kir  lifetime  analysis,  the  signal 
and  the  background  were  constrained  to  have  the  same  resolution  function;  this  allows 
the  background  to  pull  directly  on  the  fitted  signal  lifetime.  Indeed,  we  are  unable  to 
extract  unbiased  lifetimes  from  simulated  Ktt  samples  in  the  presence  of  background 
without  allowing  the  background  to  have  its  own  resolution  parameters. 


Table  3.1:  Description  of  inputs  from  the  data  and  lifetime  fit  parameters. 


Input  Data 

Description 

U 

Proper  Time 

Proper  Time  uncertainty 

PSig,i 

Signal  Probability  (determined  from  mass  fit) 

Fit  Parameter 

Description 

Signal  Lifetime 

Ssig 

Scale  factor  for  proper  time  uncertainties  (signal) 

^Sigi 

proper  time  uncertainty  for  mis-measured  signal 

/sigi 

fraction  of  mis-measured  signal  with  proper 
time  uncertainty  crsigj 

^Sig2 

proper  time  uncertainty  for  mis-measured  signal 

/sig2 

fraction  of  mis-measured  signal  with  proper 
time  uncertainty  crsigj 

TbG 

Background  Lifetime 

f TbG 

Fraction  of  background  with  lifetime 

Sbg 

Scale  factor  for  proper  time  uncertainties  (background) 

<7BGi 

proper  time  uncertainty  for  mis-measured  background 

/bGi 

fraction  of  mis-measured  background  with  proper 
time  uncertainty  ctbgi 

C^BG2 

proper  time  uncertainty  for  mis-measured  background 

/bG2 

fraction  of  mis-measured  background  with  proper 
time  uncertainty  <tbg2 

CHAPTER  4 
RESULTS 


In  this  chapter,  we  describe  the  present  the  results  of  Hkk  and  t/7r7r  measurements 
in  simulated  and  real  data.  We  show  that  based  upon  our  simulation  experiments, 
we  have  great  confidence  in  our  ability  to  extract  an  unbiased  result  from  the  data. 

4.1  Simulation 

Three  simulation  tools  are  employed  to  test  all  aspects  of  the  analysis.  We 
document  the  results  of  these  tests  here,  in  order  of  increasing  complexity. 

4.1.1  LifeFit  Internal  Generator 

In  order  to  test  LifeFit  independently  of  the  rest  of  the  analysis  chain,  LifeFit 
contains  its  own  methods  to  generate  simulated  data.  The  simulation  is  generated 
according  to  the  parameters  that  are  fit  for.  Other  fitter  inputs  such  as  and  rtii 
are  thrown  according  to  distribution  which  mock  up  what  we  see  in  the  data.  This 
allows  us  to  answer  the  question  “does  LifeFit  work?”  If  so,  LifeFit  should  be  able 
to  recover  the  input  parameters  within  statistical  errors,  and  those  errors  should  be 
correct.  The  way  we  test  this  is  to  run  thousands  of  simulation  samples  and  fits  with 
statistics  similar  to  those  found  in  the  data,  and  with  inputs  that  look  similar  to  the 
data. 

What  we  find  is  that  LifeFit  is  able  to  recover  the  input  values  of  the  toy  simu- 
lation. Figure  4.1  shows  the  pull  distribution  for  the  all  important  rsig  fit  parameter 
for  one  such  run  of  simulation.  If  we  understand  the  unbinned  fit,  we  should  see  that 
Figure  4.1  has  a width  consistent  with  one  and  a mean  consistent  with  zero.  This  is 
indeed  what  we  see. 
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Signal  Lifetime  Pull  (Internal  Generator) 

200  I 


Pull 


Figure  4.1:  Fitted  rsig  Pull  distribution  from  many  samples  of  mock  data 

internally  generated  by  LifeFit,  and  a single  gaussian  fit  overlaid. 
The  mean  of  the  fit  is  0.007  ± 0.016,  and  the  a is  1.011  ± 0.011. 
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4.1.2  FASTMC 

FASTMC  is  a derivative  of  MCFAST[26].  FASTMC  is  designed  to  be  a tool  to 
understand  charged  particle  tracking  performance.  It  allows  the  generation  of  single 
particles  and  their  decay  via  QQ.  It  also  traces  the  daughter  trajectories  through  a 
detector  described  in  an  ASCII  file,  includes  gaussian  multiple  scattering  in  the  tracing 
process,  creates  and  smears  hits  left  in  the  detector,  and  provides  a full  Kalman  fit 
to  create  reconstructed  track  candidates. 

In  this  analysis,  we  use  FASTMC  as  an  independent  infrastructure  to  test  the 
entire  analysis  chain,  including  D reconstruction  code,  in  an  idealized  environment. 
This  is  done  by  embedding  the  actual  analysis  code,  along  with  the  necessary  filling  of 
CLEO  analysis  data  structures,  inside  the  FASTMC  executable.  Since  the  unbinned 
lifetime  fitter  was  already  independently  tested  in  4.1.1,  we  can  use  FASTMC  as  a 
direct  test  of  the  analysis  code.  It  should  be  noted  that  the  detector  description  used 
was  certainly  less  than  perfect,  but  that  isn’t  important  here;  we  are  simply  testing 
the  math  in  our  analysis  code.  If  the  analysis  works,  we  should  also  find  that  a fit 
to  the  proper  time  distribution  using  our  unbinned  procedure  and  only  including  the 
signal  lifetime  and  scale  factor  (as  the  simulation  dictates)  should  result  in  a rsig 
consistent  with  the  input  value,  and  a 5sig  consistent  with  one. 

Using  a large  sample  of  ^ 1.2  million  Kit  signals  generated  with  an  input 
lifetime  of  414.95  fs  (same  as  in  QQ’s  decay.dec)  inside  FASTMC,  we  recover  a Tsig 
of  414.96  ± 0.60  fs.  We  also  find  that  the  proper  time  uncertainties  are  computed 
correctly  (5sig  of  1.0039  ± 0.0009).  This  implies  that  the  math  in  FASTMC  and  the 
analysis  code  is  very  close  to  perfectly  correct,  and  that  any  imperfections  negligible. 

4.1.3  CLEOG 

Armed  with  the  understanding  gained  from  the  Lifefit  and  FASTMC  simula- 
tions, we  use  large  samples  of  generic  GEANT  [22]-based  CLEOG  [27]  monte  carlo  to 
study  the  combined  eflPects  of  a full  detector  simulation,  backgrounds,  and  exposure 
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Table  4.1:  Fit  results  to  invariant  mass  distributions  from  CLEO  II. V monte 

carlo  for  KK,  tttt,  and  K^a.  Ktt  was  fit  to  two  gaussians  with  the 
same  mean  and  a polynomial  of  order  1,  while  KK  and  tttt  were 
fit  to  a single  gaussian  and  polynomial  of  order  1. 


Parameter 

KK 

TTTT 

Kn 

Area 

20294T177 

7668±98 

121482T381 

moiMeV) 

1864.33T0.05 

1864.50±0.09 

1864.45T0.02 

ai{MeV) 

5.26T0.05 

6.96T0.08 

5.43T0.05 

^2/(^1  + ^2) 

n/a 

n/a 

0.26±0.02 

n/a 

n/a 

1.87±0.03 

XVDOF 

175/75 

99/75 

75/73 

to  CLEO  software  on  the  measured  lifetimes  of  our  neutral  D mesons.  CLEOG  is  the 
CLEO  simulation  program  which  implements  GEANT-based  physics  and  full  detec- 
tor response  simulations  - it  is  our  best  simulation  technology.  By  the  term  “generic” 
monte  carlo,  we  mean  that  the  monte  carlo  was  generated  without  pre-selecting  cer- 
tain particle  species  from  the  truth  table.  It  is  the  most  general  style  of  monte  carlo, 
as  it  includes  all  known  backgrounds. 

The  CLEOG  simulated  data  used  in  this  analysis  is  continuum  (e+e"  — >•  udsc) 
monte  carlo,  sampling  the  available  CLEO  II.V  run  range,  and  include  backgrounds 
- this  was  achieved  by  embedding  an  selection  function  in  CLEOG  which  mocked  our 
analysis  code  and  passed  anything  which  could  possibly  end  up  in  our  analysis  chain. 
The  effective  sizes  of  samples  for  KK  and  tttt  are  over  8 times  larger  than  the  data, 
and  for  LfTT  is  « 6 times  larger  than  the  data. 

We  would  like  to  know  if  the  inputs  to  the  lifetime  fitter  and  variables  used  in 
our  selection  criteria  are  well  represented  in  the  simulation.  Plots  of  the  lifetime  fitter 
inputs,  m£),  t,  and  at  for  the  KK,  tttt,  and  Ktt  data  and  simulation  are  shown  in 
Figures  4.2,  4.3,  and  4.4.  We  observe  that  the  distributions  of  proper  time,  t,  and  its 
uncertainty,  at,  appear  to  be  simulated  quite  well;  however,  is  not.  The  results 
of  fits  to  the  Kt[,  KK,  and  tttt  invariant  mass  distributions  are  given  in  Table  4.1. 
The  input  mass  is  1.8646  GeV.  We  note  that  none  of  the  central  values  in  Table  4.1 
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Figure  4.2:  Comparison  oi  KK  proper  time,  proper  time  uncertainty,  and 

mass  distributions  in  data  and  monte  carlo.  Data  and  monte 
carlo  are  both  normalized  to  have  the  same  area. 
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Figure  4.3:  Comparison  of  tttt  proper  time,  proper  time  uncertainty,  and 

mass  distributions  in  data  and  monte  carlo.  Data  and  monte 
carlo  are  both  normalized  to  have  the  same  area. 
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Figure  4.4:  Comparison  of  Kir  proper  time,  proper  time  uncertainty,  and 

mass  distributions  in  data  and  monte  carlo.  Data  and  monte 
carlo  are  both  normalized  to  have  the  same  area. 


73 


are  consistent  with  the  input  at  the  la  level.  We  further  note  that  the  data  (see 
Table  4.7)  and  monte  carlo  disagree  about  the  central  value  of  the  mass  for  Ktt  and 
KK. 

By  using  truth  table  information,  we  can  identify  a subsample  of  correctly  recon- 
structed candidates.  Using  this  sample,  we  note  the  linear  reconstructed  mass- 
proper  time  correlation.  Figure  4.5  shows  the  proper  time  residual,  reconstructed 
minus  truth,  versus  the  reconstructed  mass  for  each  decay  channel.  This  figure  tells 
us  that  we  can  bias  our  lifetime  and  y measurements  by  reconstructing  a mean  D 
mass  that  is  not  the  truth.  Of  course,  we  can  bias  our  lifetime  similarly  by  getting 
the  decay  vertex  wrong,  or  by  misreconstructing  the  D momentum.  But  the  mass  is 
special,  in  that  we  can  do  something  about  it;  namely,  we  can  constrain  the  recon- 
structed mass  for  each  candidate  to  the  same  value.  In  the  monte  carlo,  the  value 
we  should  constrain  to  is  obvious;  the  input  D mass  (1.8646  GeV).  Figure  4.5  shows 
how  the  mass  constraint  removes  the  potential  lifetime  bias  due  to  the  reconstructed 
mass-proper  time  correlation.  Therefore,  we  use  a proper  time  calculated  with  a mass 
constraint  in  all  lifetime  fits.  Note  that  the  mass  constrained  proper  time  will  be  used 
only  in  the  signal  portion  of  the  likelihood  function;  it  is  incorrect  to  use  such  a proper 
time  for  a combination  which  is  a background. 

The  input  lifetime  to  the  simulation  is  415.0  fs.  We  look  for  biases  in  our 
reconstruction  technique  by  looking  at  proper  time  residual  and  pull  distributions  for 
Kn,  KK,  and  tttt  combinations  which  were  identified  from  truth  table  information  as 
being  correct  combinations.  If  no  bias  exists,  we  expect  that  the  residuals  and  pulls 
should  have  means  consistent  with  zero.  Further,  if  our  measured  proper  time  errors 
are  accurate,  we  should  find  the  pull  distributions  should  have  a width  consistent 
with  one.  We  show  the  residual  and  pull  distributions  fitted  to  two  gaussians  with 
the  same  mean  in  Figures  4.6  and  4.7,  respectively.  The  fit  results  are  summarized 
in  Table  4.2.  We  find  no  statistically  significant  biases. 
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Figure  4.5:  Comparison  of  mass  constrained  and  unconstrained  proper  time 

residuals  as  a function  of  reconstructed  mass  for  KK,  tttt,  and 
Kit  signal  from  CLEO  II. V monte  carlo. 
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Figure  4.6;  Proper  time  residuals  from  tagged  KK,  tttt,  and  Ktt  CLEO  II. V 

CLEOG  monte  carlo.  A fit  to  two  gaussians  having  the  same 
mean  is  overlaid. 
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Figure  4.7:  Proper  time  pulls  from  tagged  KK,  tttt,  and  Kn  CLEO  II. V 

CLEOG  monte  carlo.  A fit  to  two  gaussians  having  the  same 
mean  is  overlaid. 
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Table  4.3:  Fit  results  to  KK,  tttt,  and  Ktt  proper  time  distributions  in  CLEO 

II. V monte  carlo. 


Signal  Param. 

KK 

TTTT 

Ktv 

T-sig(ps) 

0.4144T0.0040 

0.4149T0.0058 

0.4160T0.0015 

‘5'sig 

1.016  (fixed) 

1.016  (fixed) 

1.016T0.005 

/sigi 

0.011  (fixed) 

0.011  (fixed) 

O.OlliO.002 

«^Sigi  (ps) 

0.583  (fixed) 

0.583  (fixed) 

0.583T0.070 

/sig2 

0.00057  (fixed) 

0.00057  (fixed) 

0.00057±0.00013 

<^sig2  (ps) 

8.0  (fixed) 

8.0  (fixed) 

8.0  (fixed) 

Background  Param. 

KK 

7T7T 

Ktt 

tbg  (ps) 

0.462T0.008 

0.466T0.067 

0.379T0.017 

f TBG 

0.871T0.011 

0.373T0.047 

0.894T0.028 

Sbg 

l.OOOiO.015 

1.230T0.050 

1.099T0.042 

/bGi 

0.029T0.003 

0.041T0.013 

0.155±0.018 

f^BGi  (ps) 

1.87±0.17 

2.06T0.41 

1.39T0.11 

/bG2 

0.005T0.001 

0.008T0.004 

0.014T0.003 

(7bg2  (ps) 

8.0  (fixed) 

8.0  (fixed) 

8.0  (fixed) 

Table  4.4:  Correlation  matrix  from  to  KK  proper  time  distribution  from 

CLEO  II. V monte  carlo. 


Par 

Global 

^Sig 

Tbg 

f Tbg 

‘S'bg 

/bGi 

<^BGi 

/bg2 

■^Sig 

0.31998 

1.000 

-0.238 

-0.018 

-0.003 

-0.005 

0.000 

0.010 

Tbg 

0.67941 

-0.238 

1.000 

-0.610 

-0.305 

-0.059 

0.000 

-0.103 

fTBG 

0.64233 

-0.018 

-0.610 

1.000 

0.296 

0.060 

0.000 

0.046 

Sbg 

0.61424 

-0.003 

-0.305 

0.296 

1.000 

0.524 

-0.006 

-0.103 

/bgi 

0.53683 

-0.005 

-0.059 

0.060 

0.524 

1.000 

-0.006 

-0.055 

f^BGi 

0.00693 

0.000 

0.000 

0.000 

-0.006 

-0.006 

1.000 

0.000 

/bG2 

0.17740 

0.010 

-0.103 

0.046 

-0.103 

-0.055 

0.000 

1.000 

Table  4.5:  Correlation  matrix  from  fit  to  tttt  proper  time  distribution  from 

CLEO  II. V monte  carlo. 


Par 

Global 

"^Sig 

Tbg 

fTBG 

‘S'bg 

/bGi 

(^BGi 

/bG2 

"^Sig 

0.27185 

1.000 

-0.174 

-0.019 

-0.007 

-0.004 

-0.029 

0.011 

Tbg 

0.73545 

-0.174 

1.000 

-0.651 

0.050 

-0.357 

0.031 

-0.007 

fTBG 

0.68035 

-0.019 

-0.651 

1.000 

-0.130 

0.179 

0.012 

-0.013 

Sbg 

0.37953 

-0.007 

0.050 

-0.130 

1.000 

-0.328 

0.188 

-0.038 

/bgi 

0.62644 

-0.004 

-0.357 

0.179 

-0.328 

1.000 

-0.385 

-0.066 

C’'BGi 

0.57715 

-0.029 

0.031 

0.012 

0.188 

-0.385 

1.000 

-0.380 

/bG2 

0.45624 

0.011 

-0.007 

-0.013 

-0.038 

-0.066 

-0.380 

1.000 
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To  make  sure  we  can  recover  the  input  lifetime  in  the  presence  of  background, 
we  fit  the  reconstructed  proper  time  distributions  for  each  decay  channel  with  the 
unbinned  lifetime  fitter.  Noting  that  the  three  decay  modes  in  question  are  very 
similar  kinematically,  we  can  use  the  signal  resolution  function  determined  from  the 
high  statistics  Ktt  sample  and  apply  it  to  the  smaller  KK  and  tttt  samples.  As 
noted  above,  when  evaluating  the  signal  portion  of  the  likelihood  function,  we  use 
a proper  time  and  uncertainty  calculated  with  the  neutral  D meson  candidate  mass 
constrained  to  1.8646  GeV . The  results  of  the  fits  are  shown  in  Table  4.3,  and  the 
covariance  matrices  are  given  in  Tables  4.4,  4.5,  and  4.6.  We  find  that  the  lifetime 
fitter  performs  very  well,  and  no  statistically  significant  biases  are  observed.  From 
these  results,  we  can  extract  measurement  of  ycp  from  the  monte  carlo  - it  should  be 
consistent  with  zero,  as  the  simulation  has  no  knowledge  of  mixing  in  the  D system. 
We  find 


Ukkmc  ~ 0-004  ± 0.010 

2/7T.MC  = 0-003  ± 0.014  (4.1) 

We  can  also  combine  the  monte  carlo  results  to  obtain 

Vmc  = 0-003  ± 0.009  (4.2) 

4.2  Data 

The  Ktt,  KK,  and  tttt  mass  distributions  for  all  candidates  passing  our  selection 
requirements  is  shown  in  Figure  4.8.  The  results  of  fits  to  these  mass  distributions  for 
the  signal  portion  are  given  in  Table  4.7.  Note  that  the  mean  of  the  KK  and  tttt  fits 
are  consistent  with  that  of  A'tt,  implying  that  any  problems  with  our  mass  scale  due 
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Table  4.7:  Fit  results  to  invariant  mass  distributions  from  CLEO  II. V data 

for  KK,  TTTT,  and  Kir.  Ktt  was  fit  to  two  gaussians  with  the  same 
mean  and  a polynomial  of  order  one,  while  KK  and  tttt  were  fit 
to  a single  gaussian  and  polynomial  of  order  one. 


Parameter 

KK 

TTTT 

Ktt 

Area 

2548A66 

965±37 

21517A183 

moiMeV) 

1865.09A0.15 

1864.69i0.26 

1864.98i0.05 

ai{MeV) 

5.59A0.15 

6.75i0.25 

5.65±0.14 

^2/ (^1  + ^2) 

n/a 

n/a 

0.27±0.04 

0-2 /C^l 

n/a 

n/a 

1.96i0.10 

xVdof 

96/75 

89/75 

90/73 

Table  4.8:  Fit  results  to  AM  distributions  from  CLEO  II. V data  for  KK, 

TTTT,  and  Kn.  Ktt  was  fit  to  two  gaussians  with  the  same  mean 
and  a threshold  function,  while  KK  and  tttt  were  fit  to  a single 
gaussian  and  threshold  function. 


Parameter 

KK 

'K'K 

Ktt 

Area 

2525A70 

918i36 

21286il77 

Am{MeV) 

145.432i0.007 

145.435i0.010 

145.436i0.002 

ai{keV) 

235.5i7.6 

237.6il0.3 

187.9i3.1 

^2/ (^1  + ^2) 

n/a 

n/a 

0.31i0.01 

0-2 /o-l 

n/a 

n/a 

3.35i0.10 

xVdof 

390/282 

343/282 

432/280 

to  tracking  hypothesis  {K  vs.  tt)  are  small  at  our  level  of  statistics.  For  reference, 
the  Am  distributions  are  shown  in  Figure  4.9  and  fit  results  given  in  Table  4.8. 

The  Ktt,  KK,  and  tttt  proper  time  distributions  are  shown  in  Figure  4.10. 
These  distributions  were  fit  in  an  unbinned  likelihood  fit  with  the  likelihood  function 
described  in  section  3.6.  The  fit  was  performed  as  in  the  CLEOG  monte  carlo  (see 
Section  4.1.3);  that  is,  a mass  constrained  proper  time  and  uncertainty  were  used  in 
the  signal  portion  of  the  likelihood,  and  the  KK  and  tttt  signal  resolution  parameters 
were  fixed  to  those  obtained  from  the  Kt^  fit.  The  only  difference  between  the  treat- 
ment of  the  data  and  monte  carlo  is  the  mass  to  which  the  D meson  candidates  were 
constrained. 

In  the  simulation,  the  true  mass  is  known  absolutely  {ttid  = 1.8646  GeV). 
However,  we  are  not  so  lucky  in  the  data.  According  to  reference  [29],  the  mass  of  the 
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meson  is  1864.5  ± 0.5  MeV;  i.e.,  it  is  not  very  well  known.  So  how  do  we  decide 
what  mass  to  constrain  to  in  the  data?  We  can  preserve  the  spirit  of  the  monte 
carlo  analysis  by  replacing  the  true  mass  with  our  best  estimate  of  the  mass  from  our 
data.  That  is,  we  correct  our  data  based  upon  the  simulation,  and  take  a weighted 
average  of  the  Ktt,  KK,  and  tttt  corrected  mean  masses.  The  result  we  find  for  this 
weighted  average  is  1.86514  GeV.  We  note  for  now  that  constraining  the  D meson 
candidate  mass  to  this  value  is  a source  of  possible  systematic  error  due  to  the  linear 
relationship  between  reconstructed  mass  and  proper  time  (which  have  different  slopes 
for  Ktt,  KK,  and  tttt).  The  important  issue,  which  the  mass  constraint  addresses,  is 
that  Ktt,  KK,  and  tttt  proper  times  are  all  evaluated  with  the  same  mass. 

The  results  of  the  proper  time  fits  are  shown  in  Table  4.9,  and  the  correlation 
matrices  resulting  from  these  fits  are  shown  in  Tables  4.10,  4.11,  and  4.12,  for  KK, 
TTTT,  and  Ktt,  respectively.  We  find  that  the  measured  Ktt,  KK,  and  tttt  lifetimes 
are  statistically  consistent  with  one  another,  and  that  the  resolution  functions  for  the 
signal  are  very  similar,  as  expected.  We  also  find  that  the  lifetime  of  the  background 
close  to  the  D lifetime.  From  the  correlation  matrices,  we  can  see  that  the  scale  factor 
for  the  signal  proper  time  uncertainties  and  the  lifetime  of  the  background  are  among 
the  most  correlated  with  the  signal  lifetime,  again  not  unexpected. 

From  Table  4.9,  we  can  extract  measurements  of  ycp'- 

Vkk  = -0.015  ± 0.029 

= 0.008  ± 0.043  (4.3) 

Combining  these  results,  using  the  same  prescription  used  to  combine  the  proper  time 
measurements  in  section  3.5,  we  obtain 


y = -0.008  ±0.025 


(4.4) 


Combinations/O.OOl  GeV/c 
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KK  Mass  Distribution  nn  Mass  Distribution 


1.825  1.845  1.865  1.885  1.905 

Kn  Mass  (GeV/c^) 

Figure  4.8:  Invariant  mass  distributions  from  CLEO  II.V  data  for  KK,  tttt, 

and  K-k  with  overlaid  fit  function  (see  Table  4.7  for  fit  results). 
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Figure  4.9:  Mass  difference  distributions,  m^.  — mo,  from  CLEO  II. V data 

for  KK,  TTTT,  and  Kn  with  overlaid  fit  function  (see  Table  4.8  for 
fit  results). 
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Table  4.9:  Fit  results  to  KK,  tttt,  and  Ktt  proper  time  distributions  in  CLEO 

II. V data. 


Signal  Param. 

KK 

TTTT 

Ktt 

T-Sig(p5) 

0.4106±0.0116 

0.4013±0.0169 

0.4046T0.0036 

*5sig 

1.089  (fixed) 

1.089  (fixed) 

1.089±0.015 

/sigi 

0.038  (fixed) 

0.038  (fixed) 

0.038±0.009 

^^Sigi  (ps) 

0.590  (fixed) 

0.590  (fixed) 

0.590±0.079 

/sig2 

0.00039  (fixed) 

0.00039  (fixed) 

0.00039±0.00028 

(^Sig2  (ps) 

8.0  (fixed) 

8.0  (fixed) 

8.0  (fixed) 

Background  Param. 

KK 

TTTT 

Ktt 

'Tbg(ps) 

0.436±0.020 

0.559±0.154 

0.376T0.030 

f tbg 

0.857±0.029 

0.322±0.075 

0.810T0.048 

Sbg 

1.037T0.041 

1.361±0.231 

1.171±0.075 

/bGi 

0.057T0.037 

0.358T0.123 

0.122T0.028 

C^BGi  (ps) 

0.74T0.37 

0.28±0.06 

1.34±0.16 

/bG2 

0.012±0.003 

0.042±0.013 

0.012±0.004 

<^bg2  {ps) 

8.0  (fixed) 

8.0  (fixed) 

8.0  (fixed) 

Table  4.10:  Correlation  matrix  from  fit  to  proper  time  distribution  from 

CLEO  II.V  data. 


Par 

Global 

'7'Sig 

tbg 

f TbG 

‘S'bg 

/bGi 

O^BGi 

/bG2 

^Sig 

0.30655 

1.000 

-0.214 

-0.023 

0.010 

-0.018 

-0.017 

0.006 

TbG 

0.71802 

-0.214 

1.000 

-0.625 

-0.250 

0.100 

0.222 

-0.010 

f TbG 

0.65438 

-0.023 

-0.625 

1.000 

0.259 

-0.032 

-0.099 

0.015 

Sbg 

0.48299 

0.010 

-0.250 

0.259 

1.000 

-0.175 

-0.024 

-0.008 

/bgi 

0.90744 

-0.018 

0.100 

-0.032 

-0.175 

1.000 

0.881 

0.366 

<^BGi 

0.91725 

-0.017 

0.222 

-0.099 

-0.024 

0.881 

1.000 

0.446 

/bG2 

0.47173 

0.006 

-0.010 

0.015 

-0.008 

0.366 

0.446 

1.000 

Table  4.11:  Correlation  matrix  from  fit  to 

CLEO  II.V  data. 

TTTT  proper  time  distribution  from 

Par 

Global 

^Sig 

T'BG 

f TBG 

Sbg 

/bGi 

<7bGi 

/bG2 

^Sig 

0.28551 

1.000 

-0.167 

-0.059 

0.029 

-0.011 

-0.034 

-0.009 

TbG 

0.73360 

-0.167 

1.000 

-0.633 

-0.313 

0.018 

0.233 

0.228 

f TBG 

0.67859 

-0.059 

-0.633 

1.000 

0.118 

-0.077 

-0.193 

-0.112 

Sbg 

0.73435 

0.029 

-0.313 

0.118 

1.000 

0.375 

-0.661 

0.088 

/bGi 

0.40804 

-0.011 

0.018 

-0.077 

0.375 

1.000 

-0.286 

0.061 

f^BGi 

0.67900 

-0.034 

0.233 

-0.193 

-0.661 

-0.286 

1.000 

-0.009 

/bG2 

0.30220 

-0.009 

0.228 

-0.112 

0.088 

0.061 

-0.009 

1.000 

Combinations/0.05  ps 
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KK  Proper  Time  (ps) 


7171  Proper  Time  Distribution 


-2.5  0 2.5  5 

Kti  Proper  Time  (ps) 


Figure  4.10:  Proper  time  distributions  from  CLEO  II. V data  for  KK,  tttt, 

and  Kti  with  overlaid  fit  function. 
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CHAPTER  5 

CHECKS  AND  SYSTEMATIC  UNCERTAINTIES 


In  this  chapter,  we  describe  the  various  consistency  checks  performed  to  make 
sure  our  result  is  robust,  and  provide  estimates  of  systematic  uncertainties  involved 
with  our  analysis  procedure. 

5.1  Data  Consistency  Checks 

In  our  normalized  lifetime  difference  measurements  between  KK,  tttt,  and  Ktt, 
we  expect  that  the  measurements  should  be  robust  in  the  various  ways  we  can  divide 
our  data  samples.  Put  another  way,  the  physics  of  this  analysis  does  not  change  as 
a function  of  the  kinematical  variables  of  our  D candidates,  or  when  the  data  was 
recorded,  except  due  to  statistical  fluctuations.  However,  it  is  important  to  verify 
this  by  actually  breaking  the  data  samples  up  into  pieces  as  a function  of  relavent 
variables  and  looking  for  statistically  significant  effects.  If  such  behavior  is  observed, 
it  may  be  indicative  of  detector  miscalibrations  or  anomalous  running  conditions, 
among  other  things.  Such  effects,  which  even  our  detailed  simulation  might  not  know 
about,  could  conceivably  bias  our  results.  In  the  subsections  which  follow,  we  look  for 
running  in  our  y measurements  as  a function  of  various  kinematical  variables  of  our 
D candidates  as  well  as  run  number.  If  no  statistically  significant  running  is  observed 
in  a given  variable,  we  expect  that  the  distribution  of  y measurements  will  fit  well  to 
a flat  line.  Therefore,  we  monitor  the  confidence  levels  of  fits  of  the  y measurements 
as  a function  of  different  variables.  We  observe  no  statistically  significant  running  of 
our  measurements  as  a function  of  any  of  the  variables  we  checked. 
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5.1.1  Azimuthal  Angle  (j)D 

To  see  if  we  observe  any  dependence  of  y on  azimuthal  angle  in  CLEO  II.  V data, 
we  divide  the  sample  into  four  bins  of  the  reconstructed  azimuthal  angle  of  the 
D candidate.  The  proper  time  distributions  of  these  samples  are  fit  independently 
with  the  unbinned  likelihood  fitter,  and  values  of  y extracted.  The  results  are  shown 
in  Figure  5.1.  We  do  not  observe  a statistically  significant  (f)  dependence. 
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Figure  5.1;  Mixing  parameters  and  yxK  vs.  4>d  in  CLEO  II. V data.  The 

result  of  a fit  to  these  points  to  a 0th  order  polynomial  and  its 
confidence  level  is  shown.  The  line  through  the  points  in  each 
plot  is  the  uncorrected  value  for  all  data  from  Section  4.2. 
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5.1.2  Polar  Angle  cosOd 

To  see  if  we  observe  any  dependence  of  y on  polar  angle  in  CLEO  II. V data,  we 
divide  the  sample  into  four  bins  of  cosOd,  the  cosine  of  the  reconstructed  polar  angle 
of  the  D candidate,  and  determine  independent  y measurements  from  these  samples. 
The  results  are  shown  in  Figure  5.2.  We  do  not  observe  a statistically  significant  cosd 
dependence. 
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Figure  5.2:  Mixing  parameters  y^^T^  and  yxK  vs.  cosOd  in  CLEO  II. V data. 

The  result  of  a fit  to  these  points  to  a 0th  order  polynomial  and 
its  confidence  level  is  shown.  The  line  through  the  points  in  each 
plot  is  the  uncorrected  value  for  all  data  from  Section  4.2. 
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5.1.3  Momentum 

We  check  for  momentum  dependence  by  dividing  the  entire  CLEOII.V  data 
sample  into  4 roughly  equal  pieces,  and  determine  independent  y measurements  from 
these  samples.  These  measurements  are  shown  in  Figure  5.3.  We  do  not  observe  a 
statistically  significant  dependence  on  momentum  of  the  D candidate. 
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Figure  5.3:  Mixing  parameters  y,r7r  and  y^K  vs.  po  in  CLFO  II. V data.  The 

result  of  a fit  to  these  points  to  a 0th  order  polynomial  and  its 
confidence  level  is  shown.  The  line  through  the  points  in  each 
plot  is  the  uncorrected  value  for  all  data  from  Section  4.2. 
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5.1.4  Decay  Angle  cos9* 

We  check  for  cos9*  dependence  by  dividing  the  entire  CLEOII.V  data  sample 
into  4 roughly  equal  pieces,  and  determine  independent  y measurements  from  these 
samples.  These  measurements  are  shown  in  Figure  5.4.  We  do  not  observe  a statis- 
tically significant  dependence  on  cos9*  of  the  D candidate. 


0.40 


I,  0.00 


-0.40 

0.00  0.20  0.40  0.60  0.80 

0.40 


W 0.00 


-0.40 

0.00  0.20  0.40^  0.60  0.80 

|cos(0*)| 

Figure  5.4:  Mixing  parameters  and  yKK  vs.  cos9*  in  CLEO  II. V data. 

The  result  of  a fit  to  these  points  to  a 0th  order  polynomial  and 
its  confidence  level  is  shown.  The  line  through  the  points  in  each 
plot  is  the  uncorrected  value  for  all  data  from  Section  4.2. 
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5.1.5  Vertex  Confidence  Level  prob(x^te) 

We  check  for  prob(x^jj.)  dependence  by  dividing  the  entire  CLEOII.V  data  sam- 
ple into  4 pieces,  and  determine  independent  y measurements  from  these  samples. 
These  measurements  are  shown  in  Figure  5.5.  We  do  not  observe  a statistically  sig- 
nificant dependence  on  prob(x^(3,)  ^^e  D candidate. 
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Figure  5.5;  Mixing  parameters  and  yxx  vs.  prob(x^te)  CLFO  II. V 

data.  The  result  of  a fit  to  these  points  to  a 0th  order  polynomial 
and  its  confidence  level  is  shown.  The  line  through  the  points  in 
each  plot  is  the  uncorrected  value  for  all  data  from  Section  4.2. 


94 


5.1.6  Dataset 

We  check  for  dataset  dependence  by  dividing  the  entire  CLEOII.V  data  sample 
into  6 time-ordered,  roughly  equally  sized  pieces,  and  determine  independent  y mea- 
surements from  these  samples.  These  measurements  are  shown  in  Figure  5.6.  We  do 
not  observe  a statistically  significant  dependence  on  dataset. 
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Figure  5.6:  Mixing  parameters  t/^rTr  and  yKK  vs.  integrated  luminosity; 

boundaries  between  bins  are  on  dataset  boundaries.  The  result 
of  a fit  to  these  points  to  a 0th  order  polynomial  and  its  confi- 
dence level  is  shown.  The  line  through  the  points  in  each  plot  is 
the  uncorrected  value  for  all  data  from  Section  4.2. 
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5.2  Systematic  Uncertainties 

Since  we  are  measuring  a normalized  lifetime  difference,  the  types  of  mistakes  we 
can  make  are  often  consistent  across  decay  channels  and  ultimately  cancel.  Of  course, 
this  is  not  true  in  general,  and  particularly  using  the  assumptions  we  have  made.  Some 
of  the  assumptions  we  have  made  are  that  KK,  tttt,  and  A'tt  have  identical  proper 
time  resolution  functions,  that  our  parameterization  of  the  background  is  reasonable 
and  robust,  and  that  constraining  the  mass  to  1.8646  GeV  in  the  data  is  the  right 
thing  to  do.  We  also  examine  some  equally  reasonable  reasonable  functional  forms  of 
our  likelihood,  which  amount  to  different  treatments  of  proper  time  outliers. 

5.2.1  Signal  Shape 

It  is  reasonable  to  expect  that  the  resolution  functions  for  the  signal  in  each  of 
the  decay  modes  are  very  similar  due  to  the  similar  kinematics  of  the  decays.  For 
the  background,  this  is  not  a reasonable  assumption  because  there  is  no  reason  to 
expect  the  composition  of  the  /Ctt,  KK,  and  tttt  backgrounds  to  be  the  same  (looking 
at  truth  table  information  in  the  simulation  supports  this  conclusion).  However,  the 
kinematics  of  /Ctt,  KK,  and  tttt  signals  are  not  in  fact  identical,  and  constraining 
the  shape  of  the  KK  and  tttt  signal  resolution  functions  to  that  of  i^Tr  is  a source  of 
systematic  uncertainty.  To  evaluate  this,  we  float  the  signal  resolution  parameters  of 
KK  and  tttt  (5sig,  /sigi,  crsigj,  and  /siga)  in  the  unbinned  lifetime  fit.  The  results  of 
this  exercise  form  the  data  and  monte  carlo  are  shown  in  Table  5.1. 

We  take  the  difference  from  the  high  statistics  monte  carlo  measurements  as  the 
systematic  uncertainty  due  to  our  choice  of  signal  resolution  function: 


^KK  = 0.007 
= 0.003 
5y  = 0.005 


(5.1) 
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Table  5.1:  Unbinned  fit  results  to  proper  time  distributions,  y measurements 

and  shifts  in  CLEO  II. V monte  carlo  and  data  for  Ktt,  KK,  and 
TTTT,  where  the  resolution  function  for  the  signal  in  KK  and  tttt 
were  floated. 


Sample 

'T/Cic(ps) 

'r,T7r(ps) 

T/r,r(ps) 

Monte  Carlo 

0.4171±0.0041 

0.4159±0.0060 

0.4160±0.0015 

Data 

0.4113±0.0117 

0.3975±0.0178 

0.4046±0.0036 

Vkk 

2/7777 

y 

Monte  Carlo 

-0.003±0.010 

0.000±0.015 

-0.002±0.009 

Data 

-0.016±0.029 

0.018T0.046 

-0.007±0.025 

^Vkk 

^2/7777 

Ay 

Monte  Carlo 

0.007 

0.003 

0.005 

Data 

0.001 

-0.010 

-0.001 

5.2.2  Background  Shape 

Expanding  and  our  mass  window  selection  criteria  allows  us  to  observe  how 
changing  the  composition  of  the  background  impacts  our  result.  Our  treatment  of 
the  background,  that  it  is  represented  by  two  components  - one  with  lifetime  and  one 
without,  is  put  to  the  test,  as  is  our  constraint  that  the  two  components  have  the 
same  proper  time  resolution  function.  The  results  of  varying  the  mass  window  size  are 
summarized  in  Tables  5.2  and  5.4  for  the  data  and  monte  carlo,  respectively.  Similarly, 
we  can  vary  the  overall  level  of  the  background  by  varying  the  Am  requirement; 
Tables  5.3  and  5.5  show  the  effects  of  varying  the  background  levels  by  ±50%  in  the 
data  and  monte  carlo.  We  know  our  parameterization  of  the  background  leaves  much 
to  be  desired,  and  expect  a significant  systematic  error  due  to  our  modeling  of  the 
backgound  shape. 

We  assign  the  systematic  error  from  the  monte  carlo  samples  to  avoid  sta- 
tistical fluctuations  in  the  data;  the  differences  between  the  monte  carlo  and  data 
backgrounds  are  likely  small  compared  to  the  systematic  induced  by  our  modelling. 
Varying  the  mass  and  Am  windows  are  not  the  same  systematic;  they  let  in  different 
kinds  of  backgrounds.  So  to  assign  the  systematic  error,  we  add  the  largest  differences 
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Table  5.2:  Unbinned  fit  results  to  proper  time  distributions,  y measurements, 

and  shifts  from  CLEO  II. V data  for  KK,  tttt,  and  Ktt  using  dif- 
ferent sized  mass  windows. 


Mass  Window  (MeV) 

tkk{ps) 

^TTTT  (P^) 

TjvTr(ps) 

±80 

0.4078±0.0112 

0.3984±0.0158 

0.4032±0.0034 

±60 

0.4093±0.0113 

0.3991±0.0162 

0.4040±0.0035 

±50 

0.4105±0.0114 

0.4027±0.0166 

0.4042±0.0035 

±30 

0.4138±0.0119 

0.3982±0.0174 

0.4073±0.0042 

Vkk 

VinT 

y 

±80 

-0.011±0.028 

0.012±0.040 

-0.004±0.024 

±60 

-0.013±0.029 

0.012±0.042 

-0.005±0.024 

±50 

-0.015±0.029 

0.004±0.042 

-0.010±0.024 

±30 

-0.016±0.030 

0.023±0.046 

-0.005±0.026 

^yKK 

Ay 

±80 

-0.004 

-0.004 

-0.004 

±60 

-0.002 

-0.004 

-0.003 

±50 

0.000 

0.004 

0.002 

±30 

0.001 

-0.015 

-0.003 

Table  5.3:  Unbinned  fit  results  to  proper  time  distributions,  y measurements. 


and  shifts  from  CLEO  II. V data  for  KK,  tttt. 
ferent  sized  |Am|  windows. 

and  Ktt  using  dif- 

Am 

cut  (MeV)  Tj^ft-(ps) 

^TTTT  (P^) 

'TA:7r(ps) 

±0.5 

0.4142±0.0116 

0.3899±0.0171 

0.4066±0.0037 

±1.5 

0.4052±0.0120 

0.4059±0.0170 

0.4035±0.0035 

Vkk 

J/tttt 

y 

±0.5 

-0.018±0.029 

0.043±0.047 

-0.002±0.025 

±1.5 

-0.004±0.031 

-0.006±0.043 

-0.005±0.026 

KyKK 

^UnTr 

Ay 

±0.5 

0.003 

-0.035 

-0.006 

±1.5 

-0.011 

0.014 

-0.003 

from  the  monte  carlo’s  nominal  measurement  in  quadrature,  from  both  the  mass  and 


Am  variations: 


SyKK  = 0.008 
Sy^,^  = 0.011 
5y  = 0.008 


(5.2) 
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Table  5.4:  Unbinned  fit  results  to  proper  time  distributions,  y measurements, 

and  shifts  from  CLEO  II. V monte  carlo  for  KK^  tttt,  and  Ktt  using 
different  sized  mass  windows. 


Mass  Window  (MeV) 

^/^^(ps) 

'r7r7r(ps) 

'Tirw(ps) 

±80 

0.4153±0.0039 

0.4120±0.0055 

0.4160±0.0014 

±60 

0.4148±0.0039 

0.4142±0.0056 

0.4161±0.0014 

±50 

0.4150±0.0040 

0.4149±0.0057 

0.4160±0.0014 

±30 

0.4144±0.0042 

0.4160±0.0061 

0.4155±0.0015 

Vkk 

2/7777 

y 

±80 

0.002±0.010 

0.010±0.014 

0.004±0.008 

±60 

0.003±0.010 

0.005±0.014 

0.004±0.008 

±50 

0.002±0.010 

0.003±0.014 

0.002±0.009 

±30 

0.003±0.011 

-0.001±0.015 

0.001±0.009 

^Vkk 

2/7777 

Ay 

±80 

0.002 

-0.007 

-0.001 

±60 

0.001 

-0.002 

-0.001 

±50 

0.002 

0.000 

0.001 

±30 

0.001 

0.004 

0.002 

Table  5.5:  Unbinned  fit  results  to  proper  time  distributions,  y measurements, 

and  shifts  from  CLEO  II. V monte  carlo  for  KK,  tttt,  and  Kir  using 
different  sized  |Am|  windows. 


Am  cut  (MeV) 

rKKips) 

'T,r7r(ps) 

^/^^(ps) 

±0.5 

0.4125±0.0039 

0.4178±0.0059 

0.4155±0.0015 

±1.5 

0.4111±0.0042 

0.4116±0.0057 

0.4159±0.0014 

Vkk 

2/7777 

y 

±0.5 

0.007±0.010 

-0.006±0.014 

0.003±0.009 

±1.5 

0.012±0.011 

0.010±0.014 

0.011±0.009 

^Vkk 

^2/7777 

Ay 

±0.5 

-0.003 

0.009 

0.000 

±1.5 

-0.008 

-0.007 

-0.008 

5.2.3  Treatment  of  Proper  Time  Outliers 

While  fixing  the  wider  gaussians  describing  the  poorest  resolution  signal  and 
background  to  8ps  works  well  in  the  simulation,  it  is,  after  all,  an  arbitrary  number. 
Here,  we  vary  this  assignment  to  see  what  effect  our  assumption  has  on  the  measure- 
ment oi  ycp-  We  also  dropped  the  “wide”  gaussian  altogether  and  replaced  it  with  a 
cut  on  proper  time.  The  results  are  displayed  in  Tables  5.6  and  5.7  for  data  and 
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Table  5.6;  Fit  results  to  proper  time  distributions  from  CLEO  II. V data  for 

KK,  TTTT,  and  K-k  using  different  treatments  for  proper  time  out- 
liers. 


Outlier  Treatment 

tkk{ps) 

"^TTTT  (ps) 

'T/r,r(ps) 

C’"Sig2  — <^BG2 

= 4ps 

0.4109±0.0116 

0.4012±0.0168 

0.4044T0.0036 

C^Sig2  = ^BG2 

= 12ps 

0.4108T0.0116 

0.4010T0.0168 

0.4044±0.0036 

l^il  < lOps 

0.4106±0.0115 

0.4007±0.0168 

0.4037T0.0036 

Vkk 

V'n'K 

y 

<7Sig2  = <^BG2 

= 4ps 

-0.016T0.029 

0.008±0.043 

-0.009±0.025 

C^Sig2  = <^BG2 

= 12ps 

-0.016T0.029 

0.008T0.043 

-0.008±0.025 

\ti\  < lOps 

-0.017T0.029 

0.007±0.043 

-0.010T0.025 

^Vkk 

Ay 

C^Sig2  = "^BG2 

= 4ps 

0.001 

0.000 

0.001 

(7Sig2  = <^BG2 

= 12ps 

0.001 

0.000 

0.000 

|tjl  < lOps 

0.002 

0.001 

0.002 

Table  5.7:  Fit  results  to  proper  time  distributions  from  CLEO  II. V monte 

carlo  for  KK,  tttt,  and  Kn  using  different  treatments  for  proper 
time  outliers. 


Outlier  Treatment 

'TKir(ps) 

^TTTT  (P^) 

ric,r(ps) 

<7Sig2  — C’’BG2 

= 4ps 

0.4142T0.0040 

0.4146T0.0058 

0.4158T0.0015 

0'Sig2  = CTbg2 

= 12ps 

0.4142T0.0040 

0.4148±0.0058 

0.4156T0.0015 

|tj|  < lOps 

0.4138±0.0040 

0.4145±0.0058 

0.4152±0.0015 

Vkk 

V'K'K 

V 

*^Sig2  — ^BG2 

= 4ps 

0.004T0.010 

0.003±0.014 

0.004±0.009 

t^Sig2  = ^BG2 

= 12ps 

0.003T0.010 

0.002±0.014 

0.003±0.009 

|tj|  < lOps 

0.003T0.010 

0.002±0.014 

0.003±0.009 

^Vkk 

Ay 

C^Sig2  = '^BG2 

= 4ps 

0.000 

0.000 

-0.001 

CfSig2  = ^BG2 

= 12ps 

0.001 

0.001 

0.000 

|tj|  < lOps 

0.001 

0.001 

0.000 

monte  carlo,  respectively;  we  take  the  maximum  deflection  from  nominal  in  the  data 
as  our  estimate  of  systematic  error: 


^Vkk  — 0.002 
= 0.001 
6y  = 0.002 


(5.3) 
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Table  5.8:  Fit  results  to  proper  time  distributions  from  CLEO  II. V data  and 

monte  carlo  for  KK,  tttt,  and  K'k,  where  the  proper  time  and 
proper  time  uncertainties  were  calculated  without  using  a mass 
constraint. 


Sample 

^/^/^(ps) 

^TTTT  (P^) 

'TA:,r(ps) 

Monte  Carlo 

0.4116±0.0041 

0.4153T0.0059 

0.4156T0.0015 

Data 

0.4095±0.0117 

0.4003T0.0169 

0.4054±0.0036 

Vkk 

V'K'K 

y 

Monte  Carlo 

O.OlOiO.Oll 

0.001±0.015 

0.007±0.009 

Data 

-0.010±0.030 

0.013T0.044 

-0.003±0.025 

^Vkk 

^V'K'K 

Ay 

Monte  Carlo 

-0.006 

0.002 

-0.004 

Data 

-0.005 

-0.005 

-0.005 

5.2.4  Proper  Time/Mass  Correlation  Mass  Scale 

A mass  constraint  (mo  = 1.8646  GeV  in  MC,  rriD  = 1.86514  GeV  in  data)  has 
been  employed  in  all  the  fits  performed  in  this  document  to  account  for  the  linear 
correlation  between  proper  time  and  mass  as  well  as  any  miscalibration  of  CLEO’s 
mass  scale  (see  Figure  4.5).  In  the  monte  carlo  proper  time  fits,  evaluating  the  signal 
likelihood  using  a proper  time  and  uncertainty  calculated  with  a mass  constrained  to 
the  input  mass  is  obviously  a reasonable  thing  to  do.  However,  in  the  data,  we  don’t 
know  what  the  actual  mass  of  the  meson  is.  The  world’s  knowledge  of  the  mass 
of  the  neutral  D meson  is  not  very  good;  reference  [29]  has  ttid  = 1.8645  ± 0.0005 
GeV,  which  means  we  may  have  constrained  to  the  wrong  mass  in  the  data.  In  an 
absolute  lifetime  determination,  if  one  constrains  to  the  “wrong”  mass,  each  proper 
time  measurement  and  the  resulting  lifetime  will  be  biased.  In  a normalized  lifetime 
difference  determination  like  this  one,  bias  can  also  result  because  the  slope  of  the 
lifetime  residual  versus  reconstructed  mass  is  not  the  same  in  each  decay  channel 
(again,  see  Figure  4.5). 

To  estimate  the  systematic  uncertainty  due  to  problems  with  the  CLEO  mass 
scale  as  well  as  our  choice  oi  rriD  to  constrain  to  in  the  data,  we  perform  unbinned 
lifetime  fits  on  the  data  without  using  the  mass  constraint.  On  average,  this  changes 
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Table  5.9;  Fit  results  to  KK,  tttt,  and  Kir  proper  time  distributions  in  CLEO 

II. V data,  where  the  proper  time  and  proper  time  uncertainties 
were  calculated  without  using  a mass  constraint. 


Signal  Parameter 

KK 

7T7T 

Ktt 

T-sig(ps) 

0.4095T0.0117 

0.4003T0.0169 

0.4054±0.0036 

Ssig 

1.072  (fixed) 

1.072  (fixed) 

1.072T0.015 

/Sigi 

0.033  (fixed) 

0.033  (fixed) 

0.033T0.010 

(^sigi  (ps) 

0.615  (fixed) 

0.615  (fixed) 

0.615±0.096 

/sig2 

0.00041  (fixed) 

0.00041  (fixed) 

0.00041T0.00028 

^Sig2  (ps) 

8.0  (fixed) 

8.0  (fixed) 

8.0  (fixed) 

Background  Parameter 

KK 

TTTT 

Kn 

'Tbg(ps) 

0.440T0.020 

0.511T0.109 

0.406T0.036 

f^BG 

0.852T0.031 

0.357T0.069 

0.681T0.053 

Sbg 

1.070T0.049 

1.326±0.102 

1.428±0.093 

/bGi 

0.122T0.069 

0.015T0.020 

0.103±0.031 

•^BGi  (ps) 

0.47T0.15 

3.10±1.58 

1.36±0.19 

/bG2 

0.013T0.003 

0.035T0.015 

0.012T0.004 

(^bg2  (ps) 

8.0  (fixed) 

8.0  (fixed) 

8.0  (fixed) 

the  mass  contributing  to  the  proper  time  measurement  by  amounts  varying  from 
50  to  450  keV;  which  is  on  order  of  the  potential  mass  scale  problem  implied  by  the 
detailed  simulation  (see  Table  4.1).  The  results  of  these  fits  are  given  in  Table  5.9 
for  the  data  and  monte  carlo,  and  the  correlation  coefficient  matrices  for  the  data 
fits  are  given  in  Tables  5.10,  5.11,  and  5.12  for  KK,  tttt,  and  K'k,  respectively.  For 
the  systematic  due  to  mass  scale  problems,  we  assign  the  difference  between  the 
constrained  and  unconstrained  measurements  from  the  data: 


^Ukk  — 0.005 
= 0.005 

5y  = 0.005  (5.4) 


5.2.5  Length  Scale 

Length  scale  problems,  while  paramount  in  an  absolute  lifetime  determination, 
are  less  important  a normalized  lifetime  difference;  such  mistakes  tend  to  drop  out. 
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Table  5.10:  Correlation  matrix  from  fit  to  KK  proper  time  distribution  in 


CLEO  II. V data,  where  the  proper  time  and  proper  time  uncer- 
tainties were  calculated  without  using  a mass  constraint. 


Par 

Global 

T-Sig 

tbg 

f TBG 

*5bg 

/bGi 

•^BGi 

/bG2 

"TSig 

0.32120 

1.000 

-0.219 

-0.031 

0.000 

-0.008 

0.000 

0.016 

T'BG 

0.70823 

-0.219 

1.000 

-0.636 

-0.248 

-0.048 

0.036 

-0.121 

f^BG 

0.66668 

-0.031 

-0.636 

1.000 

0.254 

0.060 

0.000 

0.069 

Sbg 

0.58827 

0.000 

-0.248 

0.254 

1.000 

0.317 

0.449 

0.059 

/bGi 

0.90073 

-0.008 

-0.048 

0.060 

0.317 

1.000 

0.888 

0.185 

f^BGi 

0.91732 

0.000 

0.036 

0.000 

0.449 

0.888 

1.000 

0.216 

/bG2 

0.27296 

0.016 

-0.121 

0.069 

0.059 

0.185 

0.216 

1.000 

Table  5.11:  Correlation  matrix  from  fit 

to  TTTT  proper  time  distribution  ii 

CLEO  II. V data,  where  the  proper  time  and  proper  time  uncer 

tainties  were 

calculated  without  using  a mass  constraint. 

Par 

Global 

^Sig 

TbG 

f'TBG 

‘S'bg 

/bGi 

<^BGi 

/bG2 

'^Sig 

0.27358 

1.000 

-0.152 

-0.094 

-0.009 

-0.011 

0.003 

-0.001 

T"BG 

0.63193 

-0.152 

1.000 

-0.532 

-0.090 

-0.217 

0.001 

0.020 

f^BG 

0.57640 

-0.094 

-0.532 

1.000 

-0.042 

0.080 

0.006 

-0.017 

Sbg 

0.34826 

-0.009 

-0.090 

-0.042 

1.000 

-0.240 

0.092 

0.044 

/bGi 

0.67483 

-0.011 

-0.217 

0.080 

-0.240 

1.000 

0.034 

-0.584 

C’^BGi 

0.30548 

0.003 

0.001 

0.006 

0.092 

0.034 

1.000 

-0.259 

/bG2 

0.64893 

-0.001 

0.020 

-0.017 

0.044 

-0.584 

-0.259 

1.000 

Reference  [25]  found  a systematic  error  due  to  length  scale  for  the  K'k  lifetime 
measurement  of  2.0  fs,  or  about  0.5%  of  a D lifetime.  From  the  data  result,  we  then 
find  the  length  scale  systematic  for  xjkk  and  y,r7r  to  be: 


^Vkk  — 0.0001 

Ar/,T7r  = 0.0003  (5.5) 


which  is  completely  negligible. 

5.2.6  Total  Systematic  Uncertainty 

The  systematic  uncertainties  are  tabulated  in  Table  5.13.  The  dominant  sys- 
tematic uncertainties  are  from  our  modelling  of  the  background  proper  time  resolution 
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Table  5.13:  Tabulation  of  total  systematic  error  for  yKK,  Umr,  and  the  com- 


bined value,  Hcp- 


Systematic  Error  Source 

yKK 

V'n'K 

Vcp 

Signal  Shape 

0.007 

0.003 

0.005 

Background  Shape 

0.008 

0.011 

0.008 

Outlier  Treatement 

0.002 

0.001 

0.002 

Mass  Scale 

0.005 

0.005 

0.005 

Length  Scale 

0.000 

0.000 

0.000 

Subtotal 

0.012 

0.012 

0.011 

Monte  Carlo  Statistics 

0.010 

0.014 

0.009 

Total 

0.016 

0.018 

0.014 

and  from  monte  carlo  statistics.  To  calculate  the  total  systematic  uncertainty  for  this 
measurement,  we  sum  in  quadrature  each  of  the  systematic  uncertainty  components 
and  find,  including  monte  carlo  statistics: 


5vkk  = 0.016 

^y-KiT  = 0.018 

5j/  = 0.014 


(5.6) 


CHAPTER  6 
CONCLUSION 

Using  9.0  fb~^  of  CLEO  II. V data,  we  have  measured  the  D mixing  param- 
eter y in  the  CP  conserving  limit  by  comparing  the  lifetimes  of  the  decays 
Ktt,  K K,Tnr.  We  find 

yxji  = —0.019  ± 0.029(stot.)  ± 0.016(syst.) 
y^^  = 0.005  ± 0.043(stat.)  ± 0.018(syst.) 

and  combined  in  a weighted  average, 

y^p  = -0.011  ± 0.025(stot.)  ± 0.014(syst.)  (6.1) 

Hkk,  VirTTi  and  the  weighted  average,  ycp^  are  consistent  with  zero,  and  has  similar 
sensitivity  to  [8].  We  are  unable  to  remove  theoretical  models  that  weren’t  already 
in  conflict  with  [12]. 

Representing  our  result  as  a 95%  confidence  level  interval: 

— 0.084  •<  yKK  0.046 
—0.086  < y^TTT  < 0.096 
-0.067  < ycp  < 0.045 

Figure  6.1  shows  this  result  in  the  context  of  other  charm  mixing  results  from  CLEO  [12] 
FOCUS  [11],  and  E791  [13].  The  figure  shows  that  the  above  measurement  is  not  in 
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D^-D*^  Mixing  Limits 


x'  (rjj/2) 


Figure  6.1:  Summary  of  some  recent  mixing  limits  in  the  neutral  D system. 
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disagreement  with  the  other  results.  Again,  for  purpose  of  illustration  we  have  taken 
6 to  be  zero  for  Figure  6.1;  this  makes  x'  — x and  y'  = y. 

FOCUS  has  advantages  over  CLEO  II. V in  proper  time  resolution  {at  = 0.08t£)o); 
it  is  so  good  they  can  use  proper  decay  time  significance  (t/at)  as  a selection  criteria 
to  assemble  a large,  fairly  clean  sample  of  candidates  that  did  not  have  a D*  tag 
(avoids  branching  fraction  penalty).  This  makes  the  resolution  of  their  ycp  measure- 
ment better  than  that  of  CLEO  II. V.  Combining  the  FOCUS  result  with  this  one 
produces 


0.010  < ycp  < 0.038 


at  the  95%  confidence  level. 

The  FOCUS  analysis,  the  E791  analysis,  and  this  analysis  can  be  combined  into 
a single  grand  CP  eigenstate  result: 

0.008  < ycp  < 0.034 

at  the  95%  confidence  level.  Graphically,  this  result  is  shown  with  the  FOCUS,  E791, 
and  Belle  results  in  Figure  6.2. 

If  mixing  is  seen  via  some  other  means  with  a high  rate,  then  our 

result  indicates  that  mixing  must  come  from  AM;  i.e.  mixing  would  take  place 
through  virtual  intermediate  states.  Our  result  is  consistent  with  the  Standard  model 
expectation  that  y is  very  small.  In  addition,  our  result  is  consistent  with  many 
extensions  to  the  Standard  Model,  which  may  significantly  enhance  AM,  but  not 
AF,  such  as  supersymmetric  models  with  quark-squark  alignment,  models  with  four 
generations,  and  models  of  light  scalar  leptoquarks  [9]. 
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Figure  6.2:  Summary  of  recent  ycp  limits  in  the  neutral  D system,  using 

eigenstate  decays. 
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